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Frequency Conversion With Nonlinear Reactance 
Chester H. Page 


\ lossless nonlinear impedance subject to an almost periodic voltage (sum of sinusoids) 
will absorb power at some frequencies and supply power at other frequencies. Necessary 
and sufficient relations among these powers are found. It is shown that simple cubic capaci- 
tors (Qo 1) are sufficient for producing any possible conservative modulation or distortion 


process 
1. General Considerations 
Let there be a set of linearly related positive frequencies w,, the relations being expressible as 


DThise,, (1) 


with the &,, rational numbers. Some of these frequencies (base frequencies) can be inde- 
pendently chosen, and the remaining ones (derived frequencies) expressed by 
Wy SS m, Whi, 


j 
or, in matrix notation, 


M2). 
It is possible to have the w, comprise subsets such that 
0Y=M?o } 


but this represents a highly artificial case in practice. 
Let the power absorbed by the nonlinear reactor at frequency w, be /?,, and for convenience 
let A a W,. 


Now 


P,Q, + WQy 
(7, + 7,M) Q,, 


and in the case where the comprise subsets: 


SP. = DBP + PPM), 


i 


But for a lossless device, >3/’?,—0 for all choices of the independent (base) frequencies, so 


WP +pPPM?)o,,=0 all 7 (Ga) 


UY +pPP Ms? —0. all 7 (6b) 
From (6a), 537 0, so power is conserved for each subset of related frequencies, as well 
as in toto. Henceforward we shall drop the superscript (j) for convenience, rewriting (6b) as 


}J,+71,M—0. (7) 











Hence the power at each base frequency is determined by the set of powers at the derived 
frequencies. Ordinarily, the base frequencies are taken as the incommensurate frequencies of 
the sources; the derived frequencies, those of the sinks. Hence the distribution of power load 
among the sources is fixed by the power dissipated in the sinks. Expanding (7) into its com- 
ponents: 

R+ 3° 2m, =0 all }. (8) 


— \ 
t 


Expressing a// w, in terms of the base frequencies: 


w= > 51, ;wn,, (9) 


j 


so that the /,,=m,, for w, a derived frequency, and /,; 
i 0 7 x) 
0 . a 
; 
for w, a base frequency. 


Then P,,=>5-4i/;,; and (8) becomes 


’ 


0= DS Aidlijt D> Ady=D Bis, (10) 
or 
DP ilij/@ S{ Puy Dol inWon 0, all }. (11) 
t i n 


Because the subscript 7 refers to base frequencies in the above, there are as many independent 
relations (11) among the /, as there are base frequencies. For the particular case of two base 
frequencies, and the /,, restricted to integers, equations (11) become equations (24) and (25) 
given by Manley and Rowe.' The result, however, is a necessary consequence of conser- 
vation of energy in any lossless system, and includes subharmonic combination frequencies. 

For further results, it is convenient to re-express the linear dependence among frequencies 


without reference to the choice of base: 


or 
‘é,,,=0 all 7 (12a) 


with all ¢,, integers. 

If there are B base frequencies, and )) derived frequencies, then M (eq(3)) is a DB matrix 
of rational numbers. Let the least common denominator of the first row of M be e,,: of the 
second row, €:9; etc., up tO epp. Then, premultiplying (3) by the nonsingular diagonal matrix, 


€11 
0 
0) 
€pp 
vields 
£,24 £,M.o, = (2, (1 >) 
£,2,+ £,2,=0, (13a) 
or 
EQ—0, 14) 


J. M. Manley and H. E. Rowe, Some general properties of nonlinear elements. I: General energy relations, Proc. IRE 44, 904-913 (1956 
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where 


D B 
E=-D| & | & 
] 
ot « 
Q) 
B\ w 


The matrix & just defined is a minimal matrix of the type E in (12); it contains no redundant 
relations among the w. We make the notational distinction because later on we shall need 
an E containing redundancies. 

Because &, is nonsingular, (13) can be solved for 


M £,', (15) 
and (7) rewritten as 


P, =P. E; 'E,. 


Let H=},£;', a 1 XD (row) matrix. Then 


7, HE, (16a) 
}j,—HE, (1 6b) 
Y—HE, (17) 
where 
D B 
B=-1 42\iZ 


-—- 


The relation (17) (with & defined by (14)) is not only necessary for conservation of power, 
but it is also sufficient for, by (14), 


> P; =pPpQ=AELQ=—0. 


t 

Equation (17) also yields insight into the nature of the frequency mixing process, for (17) 1s 
equivalent to 

P, w; > hn€nis (18) 

n 

and defining 

Pin F=hrenw; 
we have 

P=DP (19a) 


n 


2 in=0, (19b) 


sé 


so that if P?,, is considered as the power absorbed at w, in the nth “mixing process,”’ the total 
power 7’, is the net result of a set of conservative mixing processes. In any particular process, 
power is absorbed at those frequencies for which A, and e,; are of the same sign (positive 
resistance presented to the 7th source) and power is delivered at those frequencies for which 


e,, is of the opposite sign (negative resistance presented). 
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We also note that if we consider all frequencies involving a given base, the sum of the corre- 
sponding powers vanishes. For, from (2), the derived frequencies associated with the jth base 


are 


and from (8), 
R.+S\2.m,=0, 


ed 
i 
hence 
‘ 
O=.4jwnj+ >5- Bim; jo (19e) 
i 
) ae) > \p 
I, Wy) Tt > i Fai @ai, / hj 7 - y di 
1 i 


The equivalent conditions (7) and (17) are not only necessary in the generation of modulation 
and distortion products by lossless devices, they are also sufficient. That is, for any given 
M and .4, a suitable nonlinear network can be found. To prove sufficiency, we first need 


detailed relations for the behavior of a nonlinear reactor. A capacitor will be used in the 
analysis, but an inductor would be equivalent with suitable changes of notation. 


2. The Mixing Process 
Let the voltage 


Vl 


N 
24 COS (wt+4;) (20) 


be impressed on the nonlinear capacitance 


G=CY", ('>0. 
Let 8, =,/+-6, for convenience. 
Then 
Vl \ Ml \ 
( cys c( 2 cos By ) c( > 5a,(e% +4 #1) /2 ) . 21) 
Applying the polynomial law: 
v7 
“Kyp—N ATS ‘ 8 B;_1 3.\8 ' 
A, oy NA >> te ee? é P Nol? 
s 1 
where >) indicates the sum over all sets of integers s satisfying 
. 0 DJs \ 
Expanding the binomial vields 
Q=C2-* NYS May D>) e* Big 261-8 (8, — hy)! 
Let /,=2k,—s,, so that 
. \ r’ ‘ at 1:A s +/ le l; ! >» 
Q=C2-*NYD May Dy Mii /( A yl (22 
. / . ~- e ~ ” 
the /, appearing in >} being those integers between —s, and +s, that make s,+-/, even. For 


convenience, let 
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Then 


Q=C2-*NID\(Ma¥) {1 >Se"Fi/[x,+1,]}, 
8 j i a; 
and the term in braces can be rewritten as 


Ml 
EV... ef, UD... HAs, £1", 
l ly ly i=l l lu i 


Q=C2-*N155 > Maijie?F' 4i/[s, + |]. (24) 
x l 1 


making 


Now in summing over /, any set of values /,, /, lw may be paired with the set of equal 
and opposite values. Because [x,+/,] is even, the odd part of the exponential (sin >°/,8,) will 
t 


not contribute to the sum, and we can write 


Q=C2-*N!I59 > Ma‘ [s,;+1,|—' cos < Sol (w,t+6, >, (25) 
8 l i i 
where w/+86, has been restored for 8,. 
The term of frequency «, arises from values of /;, /2,. . . / that make +o,=)3/,q;, or 
Doli +(e + Io =0. (26) 


ixk 


If a// the integral linear relations among the frequencies are given by the ‘“‘mixing processes,” 


> epi =0 (EQ=0), (27 


i 


the w, terms arise from 


l, € nis ixk 
# € ph L l, 
and 
W=C'2 VNIS “SS ‘a [s,2 b> 1! asels, 4 (€ x4 1)] ' cos (w,t +6,4 > 6,0 ) 
x uo a! j 


i 


+-agk|s,+(€,x—1)]~' cos (@,t +8, > 6p /9)) . (28) 


j 


Expanding the cosines, the expression in the braces becomes 


a.” COs w,f +O) cos (re, 4 ) ls, t € ph ; 1)] _ ls ClO ps 1)] 
+-a® SID (wet -8,) sin ( Se, 6 ) [S. +(e pet 1) 7'+ [et len t+D]~ 


Reealling the definition (23) of [s+/], the coefficients in braces can be simplified, making the 


final form of Q@ become 


Q=C2*NISS SO DS)I UN ai/[s; He) } a [(se +1) Le] 
= i 


s oo 


8 +1) COS (wt +6,) cos (>) 6¢,4;)+€,« Sin (at +6,) sin > e,8))). (29) 


Differentiating for 7=Q, and writing aj*=aj**'/a, for symmetry: 


Ys \ , ‘ ‘ ‘ ‘ ] 8 
I=C2-*NYDO>5 > ea, '{ Waji/[s; +e] }ag** "[(8. +1) +e.) 
i ® = 
; € px COS (w, t+8,) sin g,—(8, +1) sin (w,t+6,) cos gy}, (30) 
be . C yf) j- 


with »,=)>0 
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If we let 


the summation 5) with S°s,;=N is the same as >) with S)s;=N-+1 and the restriction x; > 


I i 
0. But if we let s;—0, the coefficient of sin (@,¢+6,) vanishes, and [s;+¢,,|~' vanishes unless 
€ =O; but e,,=0 eliminates the cos (w,¢+6) term. Hence s;=0 contributes nothing, and 


T=C2-*NYSO SO DSSorag | { May /[s; Le pi] } {px COS (wet +6,) Sin Y,—sy SIN (wf +O) SiN gY,?, (31) 
with the integers s; >0, 5)s;=N-+1. 
i 
Noting that the coefficient of the cosine term can be summed over s’ independently of 4, 
we define 


A, Dla; Is; tép:\, (32) 
the strength of the pth process, and 
: Ou 
> \ Vo! 4 - > 
b, 2p Sella [s;te,;|=a —s 35 
making 
J=C'2 Y NYS So wa; ; A, pk COS (wf +O,) sin Pp B, sin (w,¢+@,) cos Cp}: (34) 


Pp 


The power absorbed at frequency «, is 


P= }4C2-* Nia, > 5A,j¢ x SiN ¢,, (; 


ed 35) 
p 
, CN! ; 
and if we write 17, =,..(\1, sin ¢,, 

Pj, SH, gl (36a) 

Dp 
j= HE. (36b) 

3. Order of Processes 
In eq (32) for A,, the nonvanishing terms are those for |e,;|<s}, hence >3)e,;| <S3s;=N+1 
t l 

and we define the order of a process to be S3je,,/—1. A Q=CV* device gives rise to processes 


! 
of order N, N—2, N—4, ete, because e,; must be of the same parity as s/. 

The result of a process of given order can be produced by the superposition of processes of 
other orders, introducing new frequencies for which P,;=0. That is, high-order processes can 
be indirectly vielded by low-order devices by cascaded modulation, if suitable reactance loads 
are added to the circuit and d-e bias to obtain even-order processes from odd. By introducing 
dummy frequencies involving only reactive power, the result of anv desired set of processes 
can be achieved using only simple nonlinear devices of the form Q=CV*. A constructive 
demonstration of this will constitute our proof of sufficiency. 


4. Redundancies 


The relation (36b) represents the net effect of all mixing processes, both the minimal set and 
the redundant ones. Partitioning into 


}),— HE, (37a) 
7, HE, (37b) 
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leaves redundancies in £,, making it singular, ‘so that H cannot be found from the given #,. 
In certain special cases, the frequencies are such that there are no redundant relations, and 
E=f, so that H=H=j,£;"', if the processes of E are all of order N. Then the amplitudes 
a, can be chosen arbitrarily and the A, computed (one process for each derived frequency). 
The sin ¢g, are then known, and a sufficiently large choice of C will make the g, real. 


Because 


ce Eo E,0, T E,0,, (38) 
the phases at base frequencies can be assigned, and the remainder are determined by 
0,=—E, '(#—E,.0,)=E;'e+ Me,,. (39) 


Equation (34) vields both current components, so the sink impedances are readily computed, 
In the general case of mixed-order &, the introduction of additional derived frequencies 
with corresponding null components of WJ, can convert (17) into 


p= H'E’, (40) 


with E’ of (say) third order, but the set of frequencies will usually allow redundant processes. 
Our plan of attack is to express the given & in terms of third-order processes alone, and to 
subdivide the resulting set of frequencies into subsets not possessing third-order redundancies. 
Then for these subsets, each complete Z is also a minimal &, and the subsets of mixing processes 
can be carried out by separate cubic capacitors, using common base-frequency sources, 


The only first-order process that is possible is the identity #,;=w,, because the w were defined 
to be positive and distinet. Because a third-order device can produce only first- and third-order 
processes, the problem of enumerating producible processes among a set of frequencies is 
relatively straightforward. Although a second-order device has the same simplicity, using 
Q=CV® in the sufficiency proof automatically extends the sufficiency to (a) symmetrical 
devices and (b) devices with positive slope throughout. 


5. Minimal Third-Order Processes 


Instead of choosing the source frequencies as bases, we take suitable subharmoniecs of the 
sources, so that all derived frequencies will be integral harmonics and combinations. We start 
by expressing harmonies of w in terms of minimal third-order processes, i. e., third-order 
processes defining a set of frequencies not possessing any redundant processes of third order. 

For generating an odd harmonic, 0,w, we take a sequence of odd harmonics 0, generated by 


On11=0,/3 (41) 
if integral, and otherwise 
0,4,;=(0,+1)/2, 
using the sign that makes O,,, odd. This procedure yields a sequence that can be generated 
from w by third-order processes in one and only one way. 
For an even harmonic, we again divide by 3 if possible, and repeat until we have an even 


> 


term not containing 3 asa factor. We then start an odd series with 


0, RE: - 0,=E—0, 
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and continue down from 0, to w as before, but add zero (d-c) to the sequence. In reverse, we 


find that O, and F are generated by 


20 


O,+1 


l 
hk=0,+0,-+ (zero). 


For combination tones, Aw,-+ Bw,+ . . ., we partition the sum into a sum of three combi- 
nation tones, then repartition these in the same way, ete. A partial sum of two terms can be 
partitioned by using the artifice Aw,+ Bo,=Aw,+ Bwe,+0, as was done for even harmonics. 
The original combination tone can thus be reduced to nonredundant third-order processes that 
use harmonics of the base frequencies for starting points. These harmonics can then be 


reduced as before. 
6. Choice of Base Frequencies 


Our procedure for expressing various frequencies by minimal processes was based on those 
frequencies being integral harmonics. In terms of the original arbitrary base frequencies, we 
had 

Wi = >_ Mj jn), (2) 


where m;,; are rational fractions. For each /, let the L.C.D. of the m;; be D,;,; then 


Wei = > Pi 7,/D;;, (43) 


where r;, are integers. Let the new bases be w,,=,,/D)),; so that 
Q), RQ», (), DQ». (44) 


Some of the frequencies of wz, are new frequencies (associated with zero power) and some may 
be members of the old sets w, and w,. All remaining members of w, and w, comprise the 


components of a new wp: 
Qp= FAQs. (45) 
By the minimal process procedure of the previous section, the frequencies in w,) may be associ- 


ated in sets 
Q2 SU, OQ», (46) 


e*D/ 


where m,,; are integers, so that the corresponding £, as determined by the procedure leading 
to (13) and (14) ts 


E, t lei (47) 


A given frequency may appear in several vectors wp,. 
Let the voltages at frequencies appearing in wp, be impressed on the /th cubic capacitor, 
Q=C,V*. For each of these capacitors, the relation (16a) becomes 


Hi, Pp, 


and all components are known. Assigning arbitrary nonzero amplitudes to all voltages, the 
A,, are readily computed. The phases of the base frequencies can be assigned ; for convenience 
we let each 6,,,=0. Then (39) and (35) yield 


Op, =, 


1 l\4 pp on’ 
Biri Sin - Dili jt An. 
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Using zero-impedance generators at all frequencies involved, their interconnections with the 
/ capacitors is a trivial problem. 

It has been demonstrated that any modulation-distortion operation that satisfies the neces- 
sary conditions for lossless devices can be produced with nonlinear capacitors, and in particular, 
with cubic capacitors. The demonstration of sufficiency, however, utilized perfect generators. 
It would be more satisfying to demonstrate the operation using generators at only those 
frequencies at which the net power input is positive, and using passive linear impedances at 
all sink and null-power frequencies. This can be done for any given set of base frequencies; 
variation of frequencies may lead to variation of the strengths of the several mixing processes, 
depending on the reactance-versus-frequency behavior of the elements and the particular 
circuit configuration used, 

The discussion is simplest when there is only one base frequency; all sinks are rational 
harmonics of a given source. The procedure we have used for generating minimal processes 
was such that if any base harmonic was common to two minimal processes, all lower frequencies 
involved were also common. (Except for the d-c bias needed in some processes, which can 
be handled by connecting individual bias batteries immediately adjacent to the capacitors.) 
Thus an ordered tree of generators can be assembled, with various capacitors connected at 
appropriate points on the “trunk” or upright “limbs.”’ Equation (34) gives both components 
of current at each frequency (for arbitrary assigned voltages), so impedances can be calculated 
for replacing all generators except the power source. The impedances for the null-power dummy 
frequencies will be purely reactive; those for sink frequencies will have positive resistance. By 
associating each sink with one capacitor only, all common impedances are purely reactive. 
The necessary reactances are computed from the base of the trunk upward, leaving the ‘“termi- 
nal twigs’ for sink loads and independent reactance adjustments. Such adjustments can be 
made for any given base frequency, but cannot be guaranteed for frequency variations. 

Because the general combination tones are sums and differences of harmonics of the bases, 
the harmonic-generating trees are to be interconnected with combination-tone sinks and their 
associated capacitors. In the case of harmonics that are both sink frequencies and steps in the 
development of combination tones, extra branches can be added to the harmonic-generating 
trees, so that the harmonic sink loads do not carry the harmonic current used in generating 
combination tones. Thus a// mutual impedances can be made purely reactive, and all loads 
isolated. Thus the necessary relations among frequencies and powers are also sufficient for 
realization with cubic capacitors, at any set of incommensurate spot frequencies. 


7. Regeneration and Stability 


It is apparent from eq (19b) that in a given process power is absorbed by the capacitor at 
some frequencies and delivered at other frequencies, or that in each mixing process the capacitor 
presents positive and negative resistance components at various frequencies. At a given 
frequency, some processes may absorb energy, and some may deliver energy. At a frequency 
associated with passive impedance (no generator), any power dissipated in the impedance is 
the net power delivered at this frequency by the capacitor. The total effective resistance of 
the circuit at this frequency is zero. The impressed voltage at this frequency is zero and the 
current is uniquely determined by the processes; the circuit is stable. 

On the other hand, if for some processes the capacitor delivers power at the frequency of 
one of the generators, we call it a regenerative process; it presents negative resistance to the 
generator. If the total resistance so presented is zero or negative, the circuit is sometimes 
said to be unstable. If a circuit is unstable in this sense, i. e., if it can supply power to an 


Manley & Rowe, op. cit., p. 908 
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external generator, then there is current at this frequency even in the absence of the generator. 

Hence the unstable frequency is merely any one of the frequencies generated by the network. 

If this frequency is a rational harmonic or combination tone of the other generators, the effect 

is entirely expected. But the unstable frequency may be incommensurate with all the source 
frequencies. In this case the network is said to oscillate. This condition occurs when the set 

of frequencies involved in mixing processes is such that the number of base frequencies exceeds 

the number of independent sources. Then by (19c¢) at least one frequency involving the 

excess base is associated with negative power (unless a// such frequencies are associated with 
zero power, in which case the extra base is not needed), 

This type of oscillation can occur in systems that are ordinarily considered stable, 1. e., one 
that exhibits bounded responses to all bounded driving forces, and possesses no free oscillations. 
The presence of such oscillations is, however, dependent on threshold values of generetor 
voltages, because for sufficiently small generator voltages, only integral harmonics and combi- 


nation tones can appear. 
Wasuineton, October 3, 1956 


Unpublished work, to appear in J, Wash, Acad, Sci 
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Amplitude and Phase of the Low-Frequency 
Ground Wave Near a Coastline 


James R. Wait 


\ theoretical analysis is given for the amplitude and the phase change of the ground 
wave, originating from a distant transmitter on land, as it crosses a coastline. The land 
and sea are assumed to be smooth, and homogeneous with a sharp boundary of separation. 
\ttention is focused on the effects that take place near the coastline when it is not permis- 
sible to employ arguments based on the principle of stationary phase. A limited comparison 


is made with the recent experimental work of Pressey, Ashwell, and Fowler. 


1. Introduction 


With the increased use of low-frequency radio navigational systems, it has become im- 
portant to understand the mechanism of ground-wave propagation over inhomogeneously 
conducting terrain. Theoretical treatments of this problem usually employ an idealized two- 
part media having a sharp boundary separating homogeneous media. The extension of these 
analyses to three or more media is straightforward, although tedious. In most of the previous 
work,' the numerical values of the fields are only valid for distances from the boundary greater 
than a wavelength or so, as a consequence of the approximations employed. It is the purpose 
of the present note to investigate the field near the boundary by a refinement of the stationary 
phase method. 


2. Theoretical Development 


In an earlier paper (see footnote 1), analysis and computations were presented for the 
propagation of ground waves over a mixed path on a flat earth. The boundary separating 
the two media of contrasting conductivity was assumed to be straight and approximately nor- 
mal to the propagation path. A relation was derived for the mutual impedance Z’ between 
the terminals of the transmitting and the receiving antennas, both of which are represented 
by dipoles. The conductivity of the earth to the left of the boundary is denoted by o and that 
to the right by o,. It is convenient to let 


4’ =4+4+ 44, (1) 


where Z is the mutual impedance for the two antennas, A and B, corresponding to an earth 
that is all homogeneous with conductivity ¢ and AZ is the increment to account for the portion 
of the path to the right of the boundary with conductivity o;. From the early work it follows 
that 


h hyipw I ’ l ] 
— BY¥ot+ ko F'(ro+ Ro, n) | 145 —~ (2) 
2r(7'9 T Ro) v ' ’ T iBT T Ro) B ("ort Ro)? { 


where B=27/free-space wavelength, and h, and A, are the effective heights of the antennas at 


Z 


- . . le : . a 
A and B, w=4r1077, w=angular frequency, 7=(iuw/o)*.  F(ro+-Ro, n) is a slowly varying 
Sommerfeld attenuation function defined previously (see footnote 1), and is a function of the 
separation distance ro+ > and surface impedance 7. 


1. R. Wait, Mixed path ground wave propagation: 1. Short distances, J. Research N BS 57, 1 (1956) RP2687. References to the work of 


Millington, Clemmow, Bremmer, Furutsu, and others are given here, 
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(9,7) (9%) 


land sea 





R 
A(-f5,0) B(Ro,0) 
a 
Coastline — 
Figure | Schematic representation of a land-sea boundary with a transmitter at A and a receiver at B 


The corresponding expression for the increment AZ is an integration over the media to 


the right of the boundary, which, in figure 1, corresponds to positive values of .; 


~. 2 — 
AZ = ds me f J ; rR Fir.n)F’CR ny m ( 1+ ap) 1+ “a ) cos 6 drdy, 3 
where r and ? are the distances from the variable point P to A and B, respectively, and 6 is 
the angle between unit vectors at P in the direction of increasing r and PR. The attenuation 
function F(r,n) is of the same form as that in eq (2); however, F’(2,n,n,;) is some unknown 
slowly varving function that depends on FP, n, and {= (éuw/e,)'2|. By utilizing a stationary 
phase principle, the area integration can be reduced to a line integral from r£=0 to Ry, enabling 
a solution to be carried out, To this approximation, F’(R nm) could be replaced by Fi Pn; ) 
It was admitted that this procedure leads to a result that is not valid near the boundary. 
As a further consequence of this stationary phase argument, the reflected wave from the 
boundary ts neglected. In the present instance, 2, is considered small compared with ro. It 
then can be anticipated that the important contribution of the integration will be confined to 


a region where z and y are small compared with 7». In this region it is permissible to assume: 
Fir.n)~ F (ron) ™@ Fro +h, 0), 
F’(R an) ™ 1, (subject to o; a) 
l/r 1/ro, cos 6~ (4—R,)/R, Bir+R)~Bl(rg+-r+Rh). 


In other words, the boundary is in the Fraunhofer field of the source at A. Then it follows 


that 
Bh h h 


4r“ry 
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AZ ™(n, n) 870 F'(ry, 7) (4) 





i ipr ae f - (wr Ro) l : = 
| ¢ | PR R ( +e \dyder, (5) 


R~\(r R,)?4 y* |". 


where 


It is now noted that 


, ipR ° 
| “ dy=—"- HP (8\2—R)), (6) 


which is a consequence of Campbell and Fosters? pair No. 917, which requires 8 to have a 
finite negative imaginary part. /7 is the Hankel function of the second kind. Using this 
result, it follows that 


n)iBh hy , 


<M iBhal ¢ ars F (ran), @) 
«TT I'9 


- ( 
AZ Ui 
where 


I -5| i "He Bs Rol) de. ~ 


Integrating by parts and changing the integration variable leads to 


a : > . om "BR, 7 °@—i¢ 
I=-5— 4 Hy? (BRo)—ie- ** | H(z) "de+ | HH (z)e~ "dz (9) 


for 2, >O, and 


/ 2 te BR.) . are | me—te H(2)e “ae } (10) 
ys 8Ro 


for 2yo<0, where ¢ is any arbitrarily small real number.’ 
Now it can be readily verified that 


3) 


>: ce+t?[ 1S (2) F 1H? (z)) } =et"# H(z), (11) 


so the integrations can be carried to vield 


[=i = ((1—iBR)H2 (BR) —BRH?BRo) (12) 
for 0, and 
1=i 7 (1—i8R)H.?(—BR:) +B8RHP(—BR)] (13) 


for Ry<O. After a little substitution and cancellation, it is easy to show that 
Z’ 1 AZ 
Z Z 


l+qge *e'*|( 1a) ET (a) —al] ;? (a)| (14) 


4. N. Campbell and R. Foster, Fourier integrals (D. Van Nostrand Co., Ine., New York, N. Y., 1948 
Che finite value of « insures that the integrals converge in the Riemann sens« After carrying out the evaluation, «can be set equal to zero, 
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with a=B6R,>0, and 


Y |- qe + ¢ "(1 ; 1a) tT” (a al] (a) (15 
with a=B6l, 0) In the above 
n—™ 7 —_ a : 
q ~O.O0085 ] \ ) =. o ; (16 
” ™ 0; 


where f,,- is the frequenes in megacvcles per second; q is a measure of the conductivity contrast 
between the two media and is a positive or negative real number if displacement currents in 
the ground can be neglected. Also, as a consequence of the approximate boundary conditions, 
7 <<m=120r, and therefore |AZ/Z_ is restricted to values small compared with unity. 


3. Discussion of Results 


It is instructive, at this stage, to discuss the behavior of the ratio Z’/Z for large a@ (i. e., 
large electrical distances from the boundary). The asymptotic developments for the Hankel 


function can then be emploved, to vield 


AZ : Pa\'? . = 5 
1+] ~1+ ig ( ) [1+ —— +...| (17 
js wT Sa | 2Sa- 
when Ry, >0O, and 
ity =! 1 (sea) re | 55 -* | _ 


when Ry<0. It is most significant to note that the right-hand side of eq (17), when the square- 
bracket term is replaced by unity, corresponds to the first term of a development in higher 
powers of 4 a(m/m) obtainable from the earlier paper (see appendix). That is, the formulas 
derived by the stationary phase methods are probably only valid if a>1 corresponding to 
R, large compared to the wavelength. On the other hand, the stationary phase development 
would vield a value of unity for the right-hand side of eq (18). This would also be valid at 
large distances, in wavelengths, from the boundary. The departure from unity can be inter- 
preted as a reflected wave from the boundary that forms, with the incident wave, a standing 
wave pattern on the transmitter side of the boundary. A more detailed view of the situation 
is illustrated in figures 2 and 3, where the amplitude and phase of (1+ 42/Z) are shown plotted 
as a function of Ro/A for various values of g. It is interesting to note that very near the bound- 
ary the field has a logarithmic singularity because of the behavior of the zero-order Hankel 
function as its argument approaches zero. It is believed that the solution is not strictly valid 
for very small values of a because the approximate boundary condition introduced in the earlier 
paper is not applicable at or close to a discontinuity in the surface impedance.' 

The appropriate values of the amplitude and phase 1+ 4Z/Z derivable by stationary phase 
methods are also shown in figures 2 and 3 and indicated by dotted curves. There are significant 
differences between the two methods of calculations for regions within a few wavelengths of 
the boundary. However, the phase curves on the seaside of the boundary are almost identical 
to stationary-phase evaluated case, except for regions very close to the boundary. 

It is interesting to compare these results with some phase measurements made very 
recently by Pressey, Ashwell, and Fowler... They confirmed the presence of the phase recovery 
effect that is experienced by a wave passing from land to sea. Thev also observed svstematic 


variations, about the median recovery effect, whose magnitudes are of the order of 5 very 


‘ The approximate boundary conditions are not valid for /o/A less than about |»/n9| or a small fraction of a wavelength for well-conducting 


ground, 
B. G, Pressey, G. E. Ashwell, and C, 8. Fowler, Change of phase of a low frequency ground wave propagated across i coast line, Proc 


Inst. Elee. Engrs, [B] 103, 527 (1956 
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near the coast and decay to a negligible value within a few wavelengths from the boundary 
toward the sea. These variations, which are not predicted by the present theory, are of a 
form that varied from path to path and seemed to be dependent on the topography of the 
coastline. In practically all cases, however, they observed a large and rather sudden increase 
of phase lag very near the boundary, which is in qualitative agreement with figure 3. There 
was also some evidence from Pressey, Ashwell, and Fowler’s work that a standing wave does 
exist on the transmitter side of the boundary, although the data here is rather meager. 


4. Appendix 


The stationary phase approach gives (see footnote 1) 


- ) :, d \? (® F(rota, 9) F(Ro—-2x, 1) 
F'(ro+ Ro, ",™1) Fret Ro, n)— 1€@\n m)(5 ; ) " , i = dy, 
2riB/ Jo l(rot+.r)(Ryo—.z)]2 
where 
F(r,n)=1 i2p'| e~* dse~? 
and 


o=—'¥ (2): 


If ro is large compared to 2), it is permissible to use the approximation 


Fro I.) F (ro) 
~ 
(ro +2)? yp 


and the power series expansion of the function F (22,—.+, 9) to vield 


‘ O; 
— 2i( 1 ) 
I fe 3 ( \ 0 ig 'V *P1 _ oP 
f Z \7 ri 2 4 
where 
IBR, Ur an, 
P) ) ( ) 3 


BovuLper, Co.ro., October 25, 1956. 
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Vol. 58, No. 5, May 1957 Research Paper 2757 


Vessels for the Storage and Transport of Liquid Hydrogen ' 


B. W. Birmingham, E. H. Brown, C. R. Class,’ 


and A. F. Schmidt 


The design and construction of a large-scale hydrogen liquefier by the National Bureau 


of Standards required 


containers, These containers, known 


had to be capable of road and air transport, 
This article reports on large-scale 
the NBS Cryogenic Engineering Laboratory. 


discussed, 


within the limitations of the 


1. Dewar Configuration 


The optimum Dewar configuration is spherical, 
because the sphere has the smallest area-to-volume 
ratio of any shape. In _ practice, however, such 
theoretical configurations must be tempered by the 
commercial availability of variously shaped metal 
heads. The shape of each of the NBS Dewars [1] * 
was governed by different considerations. The 440- 


liter Dewar was designed around readily available 
commercial heads. The 840-liter Dewar was de- 


signed to fit a specific space and required special 
tooling for head fabrication. Figures 1 and 2 show 
Dewar cross sections; figures 3 and 4 show the com- 
pleted Dewars. 


2. Insulation Methods 


Only two methods of Dewar insulation appeared 
feasible at the time of the design: the use of powder 
insulation, such as silica aerogel, at a moderate 
vacuum, and the use of an unfilled space at a high 
vacuum. Although the moderate vacuum associated 
with powder insulation brings many advantages 
such as reducing the importance of very minute 
leaks and making the achievement of a static vacuum 
easier—the insulation value when used in moderate 
thickness does not appear to be as great as that of 
an unfilled space at a high vacuum (up to about 
0.01-4 Hg pressure) with surfaces having a high 
reflectivity (low emissivity). As one of the primary 
goals was the achievement of a low hydrogen loss 
without making the Dewar too bulky, high vacuum 
rather than evacuated powder was chosen for the 
Dewar insulation. 

With high-vacuum insulation, heat energy flows 
to the inner surfaces by three mechanisms: thermal 
radiation, residual gas conduction, and solid con- 
duction through piping and insulating supports. 

Without any shielding, thermal radiation from 
room temperature directly to the liquid-hydrogen 


container would be excessive and would make it 


difficult to develop a low-loss hydrogen Dewar. In 
order to reduce this thermal radiation, a liquid- 
nitrogen-cooled radiation shield was planned. This 


shield completely surrounds the liquid-hydrogen 
Work supported by the U. 8. Atomic Energy Commission and the U. 8, Ait 
Force 


Now with Beech Aircraft Corp 
Figures in brackets indicate the literature references at the end of this paper 


a concurrent development of large, 
as Dewar vessels after 








| 
| 


transportable liquid-hydrogen 
the inventor, James Dewar, 


which implied ruggedness and light weight. 
liquid-hydrogen 


development conducted by 
and an 840-liter Dewar are 


Dewar 
A 440-liter 


very effort has been made to design Dewars with very low hydrogen-loss rates 
transport requirement. 


container and intercepts all ambient temperature 
radiation from the outer shell. Thermal radiation 
to the liquid-hydrogen container is then reduced by 
a factor of at least 200 and becomes of the same order 
of magnitude as the residual gas conduction and the 
solid conduction. 

When considering heat transfer to the nitrogen- 
cooled shield, residual gas conduction for a vacuum 
of the order of 0.01-4 Hg is relatively unimportant 
compared with thermal radiation. A low heat leak 


to the nitrogen-cooled shield therefore requires 
highly reflecting surfaces. Throughout the NBS 


Dewar development program high insulating vac- 
uums and surfaces of high reflectivity were em- 
phasized. 


3. Materials of Construction 


Several factors must be considered in varying de- 
grees in regard to the materials of construction of 
Dewar vessels. Mechanical properties, thermal con- 
ductivity, surface reflectivity, out-gassing rates, and 
diffusion properties are some of the more important 
factors. For equipment designed for transport, 
another factor—weight—is also of considerable im- 
portance. 

Among the various mechanical properties of metals, 
brittleness or notch sensitivity is one of the most im- 
portant in eryogenic work, because many metals 
notably the ferritic steels—become very brittle just 
below room temperature. Some of the common en- 
gineering metals whose impact strengths (as measured 
by Charpy or Izod notch tests) are satisfactory down 
to liquid-hydrogen temperatures, are austenitic steels 
(such as the stainless types, 304, 316, 321, and 347), 
copper, many copper alloys, and aluminum. 

For piping, good strength and low thermal conduc- 
tivity are necessary. Of the more common materials, 
stainless steel and Monel have these characteristics. 
Of these two, Monel not only has a thermal conduc- 
tivity approximately twice that of stainless steel, but 
it is also more expensive. Thus, in the absence of any 
other compelling factors, stainless steel is ordinarily 
used for all internal piping. 

The choice of shell materials is more difficult be- 
cause a greater number of factors are involved. For 
low weight, aluminum is indicated. For high reflec- 
tivity, both copper and aluminum are good, but 
aluminum probably maintains its reflectivity better 
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Figure 3. ,40-liter liquid h ydroge n Dewar. 


A, Liquid-level gages (N») and H,): B, pressure gages (N 
and H»); C, weight-loaded relief valve; D, insulating vacuum 
gage; E, auxiliary vacuum pump connection; F, blowers: 
Gi, shielded air-cooled diffusion pump; H, mechanical vacu 
um pump; I, liquid N, fill-and-discharge valve; J, electrical 
System switch panel and vacuum gage control unit: K. 
liquid H» discharge valve; L, liquid Hy» fill valve: M. relief 
valves (Hy) and N; 
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FIGURE 4. 840-liter liquid hydroge n Dewar. 

A, J-T valve control B, H. P. Hp» refrigerator connection; C, H. P. Hoe line filter D, L 
liquid H» fill and discharge valve: H, liquid N» fill and discharge \ 

process; K, refrigerator heat exchanger w irmup connections: L 





than copper. It is possible to use stainless steei 
covered with aluminum foil to obtain a high reflec- 
tivity. For high thermal conductivity, which is de- 
sirable in some places, either copper or aluminum can 
be used. Finally, because the solubility of hydrogen 
in aluminum is much less than that in stainless 
steel [2], aluminum should be a very good material 
with respect to hydrogen diffusion. 

For the outer shell, where low weight and high re- 
flectivity are important, aluminum would appear to 
be the ideal material. For the liquid-nitrogen-cooled 
radiation shield, low weight, high reflectivity, and 
high thermal conductivity are desirable. The latter 
factor—high thermal conductivity import- 
ant because the liquid nitrogen is in contact with only 
a restricted area of the shield, which means that ther- 
mal energy reaching the remainder of the shield must 
flow to the liquid nitrogen by solid conduction. This 
combination of requirements again indicates alu- 
minum. 

For the liquid-hydrogen container, high reflec- 
tivity, low hydrogen diffusion, and, again, high ther- 
mal conductivity are the important factors. Liquid 
hvdrogen when kept in a shell of relatively low con- 
ductivity, i. e., stainless steel, supports a temperature 
gradient that is undesirable when the vessel is sealed 
to allow pressure buildup. The top laver of the liquid 
may be as much as | to 14; deg C warmer than the 
bottom laver, and the vapor pressure in the vapor 
space above the liquid is correspondingly bigh. This 
pressure buildup is undesirable in a sealed-off con- 
tainer because it causes the safety valves to relieve 
the vessel prematurely, and in a refrigerated vessel 
causes the refrigerator to cvcle more frequently. In 
stainless or other poorly conducting shells, this strati- 
fication can be overcome by thermal equilibrators 
vertical copper or aluminum sheets or rods which 
conduct heat from the top to the bottom of the liq- 
uid—but it is simpler, if possible, to accomplish this 
by use of a good thermal-conductivity material in 
the shell itself, e. g., aluminum or copper. 

Because of the avove combination of reasons it was 
decided to fabricate most of the Dewar shells from 
aluminum. The particular tvpe of aluminum chosen 
was the weldable alloy 52S, which in liquid oxygen 
equipment had been found to be even less susceptible 
than 61S to cracking in the welds at low temperatures. 
Although aluminum had been used extensively in 
liquid-oxygen plants and Dewars, as mentioned by 
Gillette [3], this appears to be the first time its use was 
considered both for liquid-hydrogen temperatures 
and for high-vacuum vessels. The solution of the one 
foreseeable problem that arose because of this choice, 
transition joints between aluminum shells and stain- 
less steel tubing, is discussed below. 


1s 


4. Insulating Supports 


Acceleration, shock, and vibration conditions can 
be quite severe in equipment subject to road and air 
transport. For instance, the following mechanical 
service conditions were imposed as 840-liter Dewar 
design objectives: 








(1) sAeceleration: 

(a) Vertical accelerations from —3 g to +5 g, in- 
cluding acceleration of gravity. 

(b) Horizontal accelerations from —4 g to +4 g. 

Accelerations are continuously applied for periods 
not exceeding 1 sec. 

(2) Shock: 

Shocks of 5 g applied in either direction along three 
mutually perpendicular axes, with each shock being 
from 10 to 15 msec duration. 

(3) Vibration: 

(a) 0.05-in. double amplitude (total excursion) 
from 5 to 10 eps. 

(b) 0.036-in. double amplitude from 10 to 75 eps. 

Vibration continuous in any one or all of three 
mutually perpendicular planes. 

Thus it is readily seen that a Dewar capable of 
road and air transport requires a sturdy suspension 
system (i. e., a system for supporting one shell from 
another). This suspension system must also be a 
good thermal insulator in order to minimize conduc- 
tion heat leak to the cold shells. From a mechanical 
viewpoint it is desirable to have the suspension 
system rigid to eliminate the possibility of internal 
resonances that could be caused by vibration during 
transport, and from an assembly viewpoint the 
suspension system should be simple. 

For the suspension system used on the NBS 
Dewars, there appeared to be a choice between 
either stainless-steel cables or rods, or miultiple- 
contact thermal insulators which have been developed 
by the NBS Cryogenic Engineering Laboratory. 
Stainless-steel cables have been employed very suc- 
cessfully, from a thermal point of view, in liquid- 
hydrogen transport Dewars. However, their use has 
resulted in two major difficulties. First, getting the 
proper distribution of tension on all support cables 
causes a very difficult assembly problem. Second, 
experience has shown that such cables may be 
stretched considerably under repeated loading and 
thermal cycling. This stretching lowers the natural 
frequency of the suspension system and makes it 
difficult for the Dewar to meet acceleration and shock 
specifications. 

Some remarks on the development of multiple- 
contact thermal insulators by the Bureau might be 
of interest. During various experiments it has been 
observed that difficulty is experienced in transferring 
heat across metal surfaces clamped together in 
supposedly good contact in a vacuum. The surfaces 
may appear quite smooth and yet, due to microscopic 
irregularities, may present only a small area of metal- 
to-metal contact; thus, the net area for solid conduc- 
tion is quite small. This area depends on load and 
metal hardness. When one multiplies this effect by 
lamination, rigid members that are also good thermal 
insulators can be readily devised. Preliminary 
experiments by the NBS Cryogenic Engineering 
Laboratory have indicated that the apparent thermal 
conductivity of a laminated stainless steel or Monel 
insulating support under high unit loads (1,000 psi) 
would be about one-fiftieth that of the metal itself [4]. 
Further experiments [5] indicate that, by proper 
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design (dusted surfaces and moderate unit loads), 
ratios as low as one two-hundredth of the base metal 
can be achieved. 

The insulating supports developed for the NBS 
Dewars were rigid laminae formed from thin stainless 
steel (0.002 to 0.005 in. thick). Washers and disks 
were used to support axial loads, and tightly wound 
coils were used to support loads that were radial (with 
respect to the coil) (see fig. 5). 


5. Joining Aluminum to Stainless Steel 


To use metallic tubing of low thermal conductivity 
such as stainless steel in Dewars of aluminum-alloy 
shell construction, it is necessary to make aluminum- 
to-stainless-steel tubing joints that will remain 
vacuum-tight over a possible range of 130° to 

420° F. Because the coefficient of thermal 
expansion of aluminum and most of the common 
aluminum alloys is approximately 50 percent greater 
than that of stainless steel, the joint must be suffi- 
ciently ductile to withstand all stresses imposed 
upon it by thermal expansion and contraction. In 
addition, the joint must be designed to accommodate 
an amount of vibration and shock loading commen- 
surate with specified service requirements. Finally, 
the simplicity and ease of actually making the joint is 
a desirable and almost essential feature, as some form 
of field repair work may prove necessary upon occa- 
sion. 

Four basic methods of joining aluminum to stain- 
less steel have been successfully used, each of which 
employs soft solder as the filler material. (For a 
survey of methods of joining aluminum to stainless 
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steel, see [6].) In all cases, the aluminum and stain- 
less-steel tubing were initially of the same nominal 
size; however, prior to the joining operation, the 
aluminum tube was swaged out at one end for a 
length of approximately | in. so that it would fit 
easily over the outside diameter of the stainless-steel 
tube, thus forming a tubular lap joint (see fig. 6). 
In an aluminum-to-stainless-steel connection such 
as is described here, it is of definite value to utilize 
the aluminum as the outer member, because in cool- 
ing from the soldering temperature to room tempera- 
ture, and eventually from room temperature to 

420° F, the joint tends to close up and become 
even tighter than when fabricated, due to the greater 
contraction taking place in the aluminum as com- 
pared with the stainless steel. 

In brief, the four joining methods are as follows: 
(1) To eliminate the usual amount of difficulty 
involved in effectively tinning aluminum, the portion 
of aluminum to be used is first copperplated in an 
approved manner, thus providing an excellent soft 
solder base. Subsequent tinning of the copperplated 
area and the stainless-steel member provides the 
basis for a satisfactory joint. (2) An alternative 
method to that described above concerns the use of 
ultrasonic soldering equipment, which removes the 
troublesome aluminum oxide layer by cavitation, 
thereby permitting the oxide-free aluminum alloy to 
be tinned directly. (3) When it is found necessary 
to make, or repair, an aluminum-to-stainless-steel 
connection without benefit of copperplating or 
ultrasonic equipment, such as in the field, friction 
tinning can be of value. In this method, molten 
solder is applied to the aluminum working surface 
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and a sharp implement is used to rub the surface free 
of oxide beneath the solder, again permitting the 
oxide-free aluminum alloy to be tinned directly. (4) 
In order to avoid the annoying and ever-present 
possibility of oxidizing the stainless steel when mak- 
ing the joint, it is often beneficial to use a copper 
transition, or intermediate member in the connec- 
tion. The copper sleeve may be joined either to the 
inner or outer surface of the stainless steel by suitable 
silver brazing technique, and the free end is tinned 
with soft solder. The joint is made after treating 
the aluminum portion as in (1), (2), or (3) above. 

A number of other methods have been tried for 
this type of work, some of which have been erratic 
in reproducibility and others of which have met with 
little or no success. 

Possibly one of the more desirable methods of 
making a bimetallic joint involves the molecular 
bonding of aluminum and its alloys to steel. In this 
process, the aluminum is cast by any of the standard 
casting methods, against a specially prepared, clean, 
ferrous surface. In this manner, any practicable 
aluminum form can be incorporated into the joint 
design. Tests on a preliminary batch of 5 such 
bimetallic joints indicated that 4 remained vacuum- 
tight between room temperature and —320° F, and 
| leaked at any temperature. A very poor inter- 
metallic bond was observed on the one failure that 
occurred, but the remaining four appeared good. 

It appears that the best method of joining alumi- 
num and stainless steel depends greatly upon the 
application involved as well as the availability and 
suitability of equipment to do the particular job. 
Although it is not necessary to employ an expert to 
fabricate connections of this type, a certain amount 
of skill and care is of considerable value in doing good 
work. 


6. Auxiliary Design Features 


The transfer of liquid hydrogen from these large 
vessels is usually accomplished by increasing the gas 
pressure in the free space above the liquid until it is 
sufficient to force the liquid through an insulated 
transfer line to the receiving vessel at a lower 











FILLET 
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pressure. A simple, convenient method for pressur- 
izing has been developed at the Bureau. It consists 
of a circulating pump, which withdraws cold hydro- 
gen gas from the Dewar through a copper warmup 
tube and then returns the gas to the space above the 
liquid in the Dewar. In such a system, for every 
standard cubic foot of hydrogen gas circulated, about 
2.5 whr is available to cause evaporation. This 
raises the pressure over the liquid. If pressurized 
quickly enough, only the top layer of liquid hydrogen 
will come to equilibrium with the vapor. - This 
technique makes possible the transfer of subcooled 
liquid hydrogen and results in smaller losses than 
when transferring liquid hydrogen at saturation 
temperature. 

The pump used with the 440-liter Dewar is a 
“-hp oilless diaphragm paint-sprayer compressor 
with a 5.8-cfm free gas displacement. The pump is 
small and compact enough to be permanently 
mounted on the Dewar. The time required to 
pressurize the full Dewar from 1 to 15 psig is about 
3 min. The time required to pressurize the Dewar 
when half-full or nearly empty is somewhat longer. 
Transfer rates of over 2,000 liters/hr have been 
accomplished with this pressurizing technique. 

One of the predominantly experimental features 
of the 440-liter Dewar is centered around the use of 
two liquid-hydrogen lines, one at the top, the other 
at the bottom of the hydrogen container. Prior to 
fabrication, various theories existed as to the desir- 
ability of filling the liquid-hydrogen container from 
the top rather than from the bottom, and vice versa. 
If the filling takes place from the top, a great deal 
of splashing occurs, and an additional line at the 
bottom is required to empty the Dewar. If, on the 
other hand, the filling takes place from the bottom, 
flash vapor from the transfer operation bubbles 
through the liquid as it collects in the container, 
possibly entraining liquid and carrying it out the 
vent line. However, the vessel may be emptied 
through the same line. Contrary to speculation, no 
difference was noted between top and bottom ports 
in the filling operation, so it appears advisable to 
use a single valve and line leading to the bottom of 
the hydrogen container for both filling and emptying 
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The insulating vacuum of each Dewar was pro- 
vided by an air-cooled oil-diffusion pump having a 
rated speed of 260 liters/sec, together with associated 
forepump, valves, and safety features. The 840- 
liter Dewar was required to operate at ambient 
temperatures up to 130° F, at which temperature 
considerable backstreaming of oil vapor from the 
diffusion pump was found in laboratory tests. To 
prevent the oil from condensing on the cold good- 
reflective surfaces in the Dewar and thereby raising 
the emissivities, a baffle cooled by metallic contact 
with the liquid-nitrogen-cooled shield was installed 
in the pumping line to the Dewar vacuum space. 


7. Thermal Calculations 


Preliminary heat-transfer studies were carried out 
on various support and insulation methods in order 
to obtain predicted values for the heat flow into the 
liquid-nitrogen-cooled shield and the liquid-hydrogen 
container of each Dewar. As mentioned above, the 
heat flow to each of these shells can be broken down 
into three mechanisms: (1) Thermal radiation, (2) 
residual gas conduction, and (3) heat flow through 
piping and insulating supports. 

For heat-transfer calculations the Dewars may be 
considered as two separate units. The first unit 
consists of a liquid-nitrogen container at 77° K sur- 
rounded by the outer shell at 300° K; the second 
a hydrogen container at 20° K surrounded by the 
liquid-nitrogen-cooled shield at 77° K. Each of 
these units is considered separately in calculations 
because the two important figures in the results are 
the evaporative loss rates of liquid nitrogen and 
liquid hydrogen. 


7.1. Radiation 


Radiant energy transfer may be calculated from 
the Stefan-Boltzmann equation, 


Q=ceA\(Ti—T"), 


where 
(J—watts, 
o=0.0037 * 107° w in.~* deg K 
€=eMissiv ity factor, 
7’, =temperature of inner surface in degrees K, 


T’,=temperature of outer surface in degrees K, 
area of inner surface in square inches. 


The emissivity factor used in the above equation is 
obtained from a correction formula developed by 
Christiansen [7] for concentric spheres and cylinders: 





where the subscripts 1 and 2 refer to the inner and 
outer surfaces, respectively. A value of 0.04 for é 
and ¢ was used in the radiation calculations, al- 
though at the time there were no experimental data 
for the particular alloy used to substantiate this 
value. Subsequent experiments resulted in a value 
of 0.03 for these surfaces. However, irregularities 
due to piping and supports would increase the 
effective emissivity so that 0.04 appears to be a 
reasonable value. 


7.2. Residual Gas Conduction 


The residual gas conduction was calculated from 
an equation derived by S. 8. Kistler [8]. In slightly 
modified form this equation can be written, 


: M_ T,-—T, 
0.374€ 7 : ’ 
ieee |. L4(2_1)! 
/, a P 
where 
q=Wwatts per square inch of inner surface, 


(’,=specific heat of the residual gas at constant 
volume in joules per gram per degree K, 

T,, =(7T,+ T,)/2=mean temperature of the gas 
in degrees K, 


7,=temperature of the inner surface in degrees 
K, 

7,=temperature of the outer surface in degrees 
K, 

L=perpendicular distance between inner and 


outer surface in inches, 
/,=mean free path of the gas in inches at T, 
and a pressure of | mm Hg, 


P=pressure of the gas in mm Hg, 

a=accommodation coefficient factor of the 
surfaces, 

M=molecular weight of gas. 


It should be noted that g, the rate of energy flow 
per unit of area, is used directly, because at low 
pressures the concept of specific conductivity for 
the gas is meaningless. If the residual gas 1s air 
the value of /, will be /,=6.56 10-° T,, in., whereas 
if the residual gas should be hydrogen, the value of 
/, would be /,=12.2*10°°7, in. As in the radia- 
tion calculations, a correction for the Dewar con- 
figuration must be applied. It seems reasonable to 
use a Christiansen-tvpe correction factor to define 
an effective accommodation coefficient, a. 


1); 


A, and AA, are the areas of the inner and outer sur- 
faces, respectively, and @ and @ are the accommo- 
dation coefficients of the surfaces for the residual 
gas molecules. 


1 1,A/1 
ae om 


a 
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7.3. Piping and Insulating Support Conduction 


Heat flow through the piping was computed by 
yse of the standard heat-conduction equations and 
values of thermal conductivity available in the liter- 
ature [9]. Calculations on insulating support con- 
duction were made at a time when limited data 
were available on thermal contact resistance (see 
section 4). An informal communication contained 
the only data available at that time. A figure of 
approximately 0.75 mw/cm °C was given for the 
equivalent thermal conductivity of a Monel coil. 
Because stainless steel provides better insulation 
than Monel, and because the thermal conductivity 
decreases with decreasing temperature a value of 
t=0.5 mw/em °C was arbitrarily chosen for the 
calculations. 


8. Thermal Performance 


8.1. The 440-Liter Dewar 


The average total heat leak by piping, residual 
gas and support conduction, and thermal radiation 


‘for the 440-liter Dewar, during a period of time 


exceeding 1's vears of actual Dewar operation, is 
shown in table 1. 


8.2. The 840-Liter Dewar 


A representative set of performance data for the 
840-liter Dewar, over an operational period of 6 
months, is given in table 2. 

Thus it is seen that under comparable conditions 
10-> mm Hg insulating vacuum and 50° to 60° F 
ambient temperature—the hydrogen capacity, hy- 
drogen heat leak, and nitrogen heat leak of the 840- 
liter Dewar is approximately twice that of the 440- 
liter Dewar. The daily evaporation of liquid 
hydrogen expressed as a percentage of the total 


capacity is nearly the same for both vessels. With 
similar design objectives some improvement per- 


centagewise should be expected in the performance 
of the larger vessel. In this case, however, more 


‘ stringent mechanical service conditions resulted in 


more heat flow through insulating supports and 
vullified such improvement. 


TABLE l., 


Liquid hydrogen 


Outer 


Residual shell Daily evapora- 


pressure ® | tempera Heat leak tion (percent- 
ture age of total 
capacity) 
mm He F 7 a / 
1o-* 60 3.7 2. 29 
5x10 “0 2.0 1. 24 
10 i 1.3 (1.4 bO. SL (O89 
10 mw >O.Y (OLS) b. 56 (0. 46) 


* For the purpose of studying the effect of a variable insulating vacuum, the tabulated data for pressures above 10 
> mm Hg was usually maintained during normal operation. 


ak to the system, because a vacuum between 10 and 10 
>» Numbers in parentheses represent design calculations. 








Without additional shields 


9. Refrigeration 


It is frequently desirable to store liquid hydrogen 
without loss. For short periods this may be accom- 
plished by sealing the hydrogen container and allow- 
ing a moderate pressure rise. In Dewars of the type 
discussed here, storage durations of 1 week or more 
result in a pressure increase of about 1 atm within 
the hydrogen system when the hydrogen container 
is initially about 90 percent full of liquid. 


TABLE 2, Performance data for 840-liter Dewar 


Liquid hydrogen Liquid nitrogen 


Outer shell ; 
tempera- 
ture 


Daily evap- 
oration (per- 


Residual 


pressure Daily evap- 


Heat leak centage of Heat leak oration 
total 
capacity) 
mm Hg F u w/ u Liters 
10 52 2.5 0. 80 118.0 63. 2 


To keep hydrogen without loss indefinitely, it is 
necessary to periodically refrigerate the hydrogen to 
prevent excessive pressure rise in the hydrogen shell. 
Several refrigeration schemes have been proposed 
and developed, which operate on the principle of 
recondensing the vapor formed from the liquid hydro- 
gen evaporation by circulating a colder gas, used as 
a refrigerant, through a closed condensing coil. The 
Bureau has developed a refrigerator that utilizes 
the principle of reliquefving the vapor formed from 
evaporation of liquid hydrogen by circulating the 
hydrogen vapor in a simple Joule-Thomson-type 
hvdrogen liquefier. A high refrigeration efficiency 
is achieved. The components of this refrigeration 
system include the usual compression, purification, 
purging, and control equipment, along with the 
necessary heat exchangers and expansion valve. 

The 840-liter Dewar is capable of being refrigerated 
by equipment using this new application of the 
Joule-Thomson process. The only components op- 
erating at low temperatures, and thus requiring in- 
sulation, are the heat exchangers and expansion 
valve. These were mounted in the Dewar insulating 
vacuum space. The remainder of the refrigeration 


Average total heat leak 


Liquid nitrogen 


With additional shields 


Heat leak Daily Heat Daily 
evaporation leak « evaporation 
uw Liters u Liters 
(63.8 (34. 2) 
662. 5 (61.1) b 33. 5 (32. 7) 52.0 27.8 


‘mm Hg were obtained by introducing a small 


It will be noted that the addition of two rigidly supported aluminum radiation shields [9, fig. 1] between the outer shell and the liquid-nitrogen-cooled radiation 
Shield reduced the actual heat leak to the liquid nitrogen shield approximately 16 percent 
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equipment is mounted on a dolly, which can be 
brought easily to the Dewar and connected. Con- 
struction of the refrigerator in this manner has re- 
sulted in maximum efficiency and simplicity because 
all connections are warm and the Dewar insulating 
vacuum is utilized to insulate all low-temperature 
components of the refrigeration cycle. 

A detailed description of this reliquefying hydrogen 
refrigerator is planned for a subsequent report. 


10. Conclusions 


The National Bureau of Standards has developed 
rugged, lightweight vessels for storing and transport- 
ing liquid hydrogen. These vessels have a very low 
loss rate. Techniques employed in this development 
include vacuum-tight aluminum welding, joining 
aluminum to stainless steel, the use of novel insulat- 
ing supports, and methods for improving the re- 
flectivity of metallic surfaces. It appears entirely 
feasible, by further refinements in producing surfaces 
of high reflectivity and emploving mechanical designs 
with less support conduction, to effect a still further 
improvement, making possible the low-loss trans- 
portation of liquid helium. 


11. References 


Government erects $3,500,000 low temperature laboratory, 
Refrig. Eng. p. 1082 (1953). 
[2] Donald P. Smith, Hydrogen in metals (Univ. Chicago 
Press, 1948). 

[3] H. W. Gillette, The behavior of engineering metals, p. 353 
(John Wiley & Sons, Inec., New York, N. Y., 1951). 

{4] Russell B, Scott (private communication). 

[5] R. P. Mikesell and R. B. Scott, Heat conduction through 
insulating supports in very low temperature equipment, 
J. Research NBS 57, 371 (1956) RP2726. 

[6] Morton C, Smith and David D. Rabb (private communi- 

cation). 

[7] Max M. Jakob and G., A. 
transfer, Ist ed., p. 125 
New York, N. Y., 1942). 

S. S. Kistler, The relation between heat conductivity and 
structure in Silica Aerogel, J. Phys. Chem. 39, 79 
(1935). 

[9] Robert L. Powell and William L. Blanpied, Thermal con- 

ductivity of metals and alloys at low temperatures, 

NBS Cire. 556 (1954). 


heat 
Ine., 


Hawkins, Elements of 
(John Wiley «& Sons, 


[8 


12. An Improved Vessel 
A. F. Schmidt and B. W. Birmingham 
Subsequent to the work reported in the preceding 
sections, the 440-liter liquid-hydrogen Dewar was 
converted in capacity and design to a 500-liter 
Dewar with several predominant modifications and 











improvements. A comparison of figures 1 and 7 wil] 
service to illustrate these changes. In_ brief, the 
aluminum 440-liter spherical container was replaced 
with a 500-liter stainless-steel dish-ended cylindri. 
cal tank, externally silver-spraved in order to ip. 
crease the surface reflectivity. The vertical stainless. 
steel center support tube was decreased in wall! 
thickness to 0.016 in. (one-fifth of its former value). 
and lateral motion of the inner container with re. 
spect to its radiation shield was nullified by 2 sets 
of 3 stainless-steel aircraft cables, *s2 in. in diameter. 
which replaced the stainless-steel coil strips pre- 
viously used; \-in. aircraft cables also replaced the 
entire bottom coil assembly. The rigidly supported 
double set of ‘floating radiation shields’? was ex- 
changed for a single ‘floating shield”’ attached to the 
outer shell by 3 long stainless-steel tubes, and Teflon 
bumpers spaced the shield from the nitrogen shel] 
at the top. In order to reduce the over-all height of 
the assembly, vent lines and valves were brought 
through the side of the Dewar. 

The thermal performance of this 500-liter Dewar 


over a l-month period is given in table 3. 


TABLE 3.— Thermal performance of 500-liter Dewar 


Liquid hydrogen Liquid nitrogen 


Residual Outer shell 


pressure termpera- Daily evapo- Daily 
ture Heat ration (percent Heat evapora- 
leak age of total leak tion } 
capacity 
mm He I ” o " Liters 
lo un 0.47 0), 26 29.5 15.8 


A note of interest may be inserted here relative to 
one aspect of safety involved in the use of this vessel, 
Users have been concerned about the consequences 
of sudden loss of insulating vacuum for containers} 
of this type. A dangerous rise of pressure was feared 
In this container a 1-in.-diameter hole was suddenly 
opened in the outer wall of the Dewar vacuum space, 
resulting in the entire contents—both hydrogen and 
nitrogen—evaporating and being vented through the, 
respective %-in. vent lines. The pressure rise result- 
ing from the rapid evaporation did no damage. 

In conclusion, these changes confirm that applica- 
tion of design modifications recommended in section 
10 were justified. A liquefied-gas-storage vessel has 
been developed with a daily hydrogen evaporation 
rate of approximately \ percent of the total capacity 
in the 500-liter size range. It is planned to evaluate 
this vessel as a liquid-helium-storage Dewar. 
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FIGURE 7. 


500-liter liquid hydrogen Dewar. 


Bovu.per, Coro., October 27, 1955. 
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Infrared Spectra of Polychlorobenzenes 


Earle K. Plyler, Harry C. Allen, Jr., and E. D. Tidwell 


The infrared spectra of 12 chlorinated benzenes have been measured from 
In addition, the spectra of benzene, monofluorobenzene, monobromobenzene, and 
An assignment of the fundamental bands of the 
Some of the frequencies of the other compounds 


microns. 
monoiodobenzene have been included. 
monosubstituted benzenes has been made. 
have also been assigned. 


2 to 38 


The spectrum of each compound is shown in the figures, and a table containing all the 


observed bands is included. 


1. Introduction 


There is now available the complete group of 
substituted chlorobenzenes, from monochloroben- 
zene to hexachlorobenzene. These have been puri- 
fied by the Hooker Electrochemical Co. and have 
been made available through R. H. Kimball. <A 
measurement of the infrared spectra of these chlorin- 
ated benzenes will serve two purposes. First, it 
will serve as a catalog of the spectra of these com- 
pounds, which will be useful for identification and 
for intercomparison with other substituted benzenes. 
Second, some of the typical vibrations can be traced 
from one substituted molecule to another, thus 
showing the effect of position of the chlorine atoms 
in the molecules on the vibrational frequencies. 

It would be a gigantic task to try to assign all the 
observed bands for 12 chlorinated benzenes. This 
has not been undertaken, but a study of some of the 
types of vibrations has been made and the classifica- 
tions listed. 

In addition to the spectra of the chlorobenzenes, 
there have been included the spectra of monofluoro- 
benzene, monobromobenzene, and monoiodoben- 
zene. These vibrational bands of monosubstituted 
benzenes have been classified and the wave numbers 
compared with those of monofluorobenzene. 


2. Experimental Methods and Results 


The spectra were measured on a prism spectrom- 
eter. Prisms of NaCl, KBr, and KRS-5 were used for 
the measurements from 2 to 38 uw. Cell thicknesses 
from 0.05 to 1.6 mm were used in different parts of 
the spectrum. The cell thickness was varied in 
regions of strong absorption so that the position of 
band center could be adequately determined. Where 
a thinner cell was used, the results are placed on the 
graphs as an insert curve. Several of the chlorinated 
benzenes are normally solids and have been meas- 
ured in solution. Carbon tetrachloride was used for 
the solvent from 2 to 8 u, carbon disulfide from 8 to 
22 and benzene or methylevclohexane from 22 to 


Ze ft, 
38 u. In some regions the absorption of several of 





the liquids was so intense that the spectrum was 
determined in solution. The bands measured under 
these conditions are designated on the figures. 
Carbon tetrachloride has an absorption band of 
medium intensity in the region of 6.4 uw, and bands of 
the solutes in this region could be determined only 
approximately. The absorption spectrum of the 
compounds in the CCl, solution show the CCl, band 
at 6.4 u plus any absorption of the solute. 

The observed spectra of the chlorinated benzenes 
are shown in nine figures. The absorption spectra of 
benzene, fluorobenzene, chlorobenzene, bromoben- 
zene, and iodobenzene are shown from 2 to 38 yu in 
figures 1 and 2. Figures 3 and 4 include the spectra 
of o-dichlorobenzene, m-dichlorobenzene, and p-di- 
chlorobenzene. The p-dichlorobenzene was measured 
in solution. The spectra of the three trichloro- 
benzenes are shown in figures 5 and 6, and the spectra 
of the tetrachlorobenzenes are included in figures 8 
and 9. The spectra of pentachlorobenzene and hexa- 
chlorobenzene from 2 to 15 uw are shown in figure 7 
and from 15 to 38 yu in figure 9. 

In the region from 23 to 38 u the spectrum of hexa- 
chlorobenzene is represented by a broken line. In 
this case, and other spectral regions of the com- 
pounds, the broken line indicates that the spectrum 
is less accurately determined than for the regions of 
the solid line. The uncertainty arises either from 
atmospheric absorption or some absorption present 
in the solvent. Bands of low intensity in the regions 
of solvent absorption are probably not observed. 

Table 1 is a list of the observed bands of all the 
compounds. The wavelength, wave number, and in- 
tensity are listed for the bands of each substance. 
Bands of low intensity are not included. The in- 
tensity is given in terms of percentage of transmis- 
sion, and is included primarily for identification pur- 
poses. No correction has been applied for radiation 
loss at the cell windows, nor has the percentage of 
transmission been reduced to a uniform cell thick- 
ness, but sufficient data are given so that a quantita- 
tive comparison of intensities can be made. It is 
estimated that the wavelengths are correct to +0.03 
uw in the region from 2 to 15 uw and to +0.06 uw from 
15 to 38 up. 
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Figure 1. Infrared absor plion spectia of benzene, fluorobenzene, chlorobenzene, bromobenzene, and todoben zene from 2 to 15 microns. 
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3. Discussion of Results inte iy i ee 
“— a Fi 

Any complete assignment of the vibrational fre- | HY a" 
quencies as determined from liquids or solutions is | . 
very difficult because the distinctive band shapes ~ hd ~— 
arising from the rotational fine structure are washed | a 
out. In the present case the difficulty is aggravated | . ue 
because the spectra beyond 38 yu were not obtained. | ; , 
Thus many of the important low-lying fundamentals —_ 


-_ 


were not observed. However, with the help of the 
excellent work of Randle and Whiffen,’ some of the | 
bands have been assigned. Some help was also | 083 to ante 
obtained from the calculations of Garg.” , ' 





R. R. Randle and D. H. Whiffen, Molecular Spectroscopy Rept. of Conf 
Hydrocarbon Research Group, Inst. Petroleum (Oct. 1954 
S. N. Garg, Current Sci. 22, 208 (1953 


Figure 2. Continuation of the spectra of the compounds of 
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Infrared Spectra of Chromatographically 
Fractionated Asphalts 


James E. Stewart 


Undegraded and weathered specimens of three roofing asphalts were chromatograph- 


ically 
absorption spectra of the fractions, 


separated into three components after removal of the asphaltenes. 
recorded in the 


The infrared 
2- to 15-4 spectral region, provide infor- 


mation about the chemical constitution of the material, the differences among asphalts from 
different sourees, and the changes in constitution produced by weathering. 


1. Introduction 


Work on the constitution of asphalt and the de- 
gradation of asphalt on weathering has been under 
the National Bureau of Standards for some 
time. This work has been facilitated greatly by the 
development) of chromatographic techniques [1] ? 
that permit the separation of structurally similar 
components for chemical and physical study. One 
phase of this study involved the use of infrared 
spectroscopy in an attempt to elucidate the chemical 
structure of the asphalt components and to deter- 
mine changes that occur in the components when the 
asphalts are weathered Infrared spectroscopy [2,:3] 
has been useful in studving the chemical steucture of 
chromatographically fractionated petroleums of the 
tvpe used as extender oils in synthetic rubber. It 
seemed reasonable, therefore, to apply infrared 
methods to the study of asphalts as well. 

The purpose of this report is to present the results 
of a brief study of three airblown roofing asphalts of 
different origin. These asphalts are product 3. of 
asphalts 1, II, and IIL, the properties of which were 
reported in detail by Greenfeld [4]. The asphalts 
and their components will be referred to as 1, IIT, and 
II], respectively, in the following discussion. Both 
undegraded materials and specimens degraded by 30 
months of natural weathering were examined. 


wil at 


2. Experimental Procedures 


The asphalts were separated into four groups of 
components by L. R. Klemsehmidt of the NBS 
Floor, Roof, and Wall Covering Section, Normal 
pentane was used for the primary separation of the 
asphalts into asphaltenes and maltenes. The mal- 
tenes were fractionated further [1] into three groups 
of components, designated as water-white oils, dark 
oils, and asphaltic resins, by eluting the n-pentane 
soluble portions of the asphalts from a fuller’s-earth 
column, with nv-pentane, methylene chloride, and 
methyl ethyl ketone 

The distribution of components in the’ three 
asphalts, both before and after 830 months of exposure 


outdoors at the National Bureau of Standards, is 
This work was supported ntly by the Asphalt Roofing Industry Bureau 
1 National Bureau of Standard 
Figures in bracket licate the literature references at the end of this pap 


given in table 1. The distribution of components 
changes during weathering through the conversion of 
certain materials into water-soluble and_ volatile 
degradation products and by the alteration of ad- 
sorption properties because of chemical reactions. 

Infrared absorption spectra were obtained with 
® Perkin-Elmer model 21 spectrophotometer, equip- 
ped with a sodium chloride prism. In order to 
avoid the interference of solvents it was decided to 
observe the material as liquid films. A cell was 
constructed that could be heated to about 60° C 
with a nichrome heater to facilitate loading the vis- 
asphaltic-resin fraction. Close control of 
temperature was not required because of the in- 
sensitivity of the absorption by these materials to 
temperature changes. The windows were separated 
by a lead spacer 0.05 mm thick, cut away at the top 
and bottom. Aslot was cut in the top of one window. 
Samples were introduced into the slot and permitted 
to run into the cell, sometimes with the aid of suction 
applied at the bottom of the cell. After a spectrum 
was recorded, the sample was removed by suction 
and the cell rinsed with solvent. To assure con- 
stancy of the cell thickness, it was not disassembled 
during the series of measurements. 

In addition, the spectra of the asphaltenes were 
obtained by suspending the material in potassium 


COUSs 


; ; - : ‘ om 
iodide pellets and in mineral oil. The use of 
TABLE 1 Analyses of asphalts 

j ! 

Pyvye I If III 

Unweathered material, stored 30 months in covered container 
Asphaltenes 41.2 41.0 39.2 
Water-white oils 27.7 25. 5 23.6 
Dark oils 7.5 20.5 21.8 
Asphaltic resins SS 10.3 10.5 

Weathered, 30 months ® 

Asphaltenes 11.6 42.0 43.6 
W ater-white oils 18.7 20.1 19. 1 
Dark oils 8.5 14.6 14.6 
Asphaltic resins v0 9.4 8.5 
Weight loss 19.7 11.7 11.5 


distribution of components based on the weights of the original 
n thickness and 12 in.? in area 


* Percentage 


coatings, which were 0.02 in 








potassium iodide provided spectra that free 
from losses by seatter, in contrast with those ob- 


were 


tained by using potassium bromide. However, 
because of adsorbed water, alkali halide pellets 
absorb in the 3u region where O--H groups are 
observed. The use of mineral-oil Suspensions per- 


mits observations in this important region 


3. Infrared Spectra 


The value of infrared spectroscopy in studies like 
the present from the fact that certain 
functional groups absorb infrared radiation in fairly 
narrow regions of wavelength, relatively independent 
of other functional groups in the molecule. Hence, 
in a mixture of many different compounds, common 
structural often be determined For 
example, all of the asphalt fractions near 
3.40, 3.49, 6.82, and 7.2 uw because of the stretching 
and angle-deformation vibrations of the CH. and 
CH, groups of alkanes. On the other hand, many 
functional groups have no characteristic absorption 
or absorb so weakly that they cannot be detected 
C Sand S—S groups are of this nature. Further- 
more, chemically different functional groups oftes 


one 


arises 


groups can 


absorb 


have similar vibations and therefore absorb at a 
common wavelength. Conjugated carbon chains 
and aromatic rings, for example, have bands at 


6.25 4. In the analvsis of the fractions that follows, 
liberal use has been made of the standard reference 
books in chemical infrared analysis; ino particular, 
the book of Bellamy [5] has been invaluable 











Kigtre ll, Jafrared spectra of u 


3.1. Water-White Oils 


The water-white oil fractions (fig. 1) are seen to 
he predominantly saturated vliphatics, with a slight 
aromatic content producing a ring vibration at 
6.25 uw. Very likely the aromatic content is largely 
of a paradisubstituted nature leading to the absorp- 
tion band at 12.2 u. From the relative intensities of 
the 6.25- and 12.2-u aromatic bands from sample 
to sample, it appears that samples Il and TIT have 
about the same aromatic content, whereas sample 
| is more aromatic by a factor of about 2 or 3.) An 
idea of the mean chain length, and thus of branching. 
is Obtained by observation of the intensities of the 
(—H band at 13.8 « and the shoulder found on the 
short wavelength side of this band. In general. the 
longer chains absorb at longer wavelengths than the 
shorter Simple II evidence of 
slightly more branching than sample I], and sample 
I is considerably more branched. This observation 
is verified by the intensities of the CH, bands at 
7.24 uw, the more highly branched samples showing 
stronger CH, absorption 

The presence of eveloparaffins in 
mixtures is believed [2] to result) in 


chains. shows 


hy drocarbon 
increased ab- 


sorption in the rather broad region centered at 
about 10.2 4. Henee, sample IIT is slightly more 
eveloparaffinic in nature than IT, and sample I is 


much more evcloparaffinic 

The water-white oil separated from a weathered 
sample IL is almost identical to the unweathered 
sample, except for a slight decrease in intensity of 
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the aromatic bands and perhaps a slight reduction 
in branching. Sample II] displays, in addition to 
reduced aromaticity, a more pronounced reduction 
in branching as well as reduced absorption in the 
eveloparaffin region After weatherme, the water- 
white oils from asphalts II] and I] are nearly ident- 
eal, with the former being very slightly less branched. 
The aromatic and evcloperaffin absorption of sample 
| is reduced after weathering, but the branching 
seems to be unaltered 

The relatively vreater loss of eveloparaffin and 
branched aliphatic components during weathering 
is in agreement with the evreater reactivity of such 
Furthermore, sample 1, which ts richest 
loses a greater fraction of the 


table | 


compounds, 
in these compounds, 
water-white oils originally present (see 


3.2. Dark Oils 


The dark oils (fig. 2) are highly aromatic im nature, 
es shown by the bands at 6.21, 11.48, 12.24, end 
13.35 uw. Substitution on the rings in various ways 
or fused rings produce the three-band pattern be- 


iween Tl and 13.5 uu. In addition, two bands are 
found at 13.68 and 13.87 u that are difficult to in 
terpret. They are probably due in part to CH, 


chains, as in the water-white oils, because there are 
ertainly aliphatic CH, and CH, groups present. 
Bands comparable in appearance to these were not 
found in the heavy petroleum fractions previously 


FiGuRE 2 Infrared 


i, Unweather 


I] 





d 


studied [2, 3}. 
the form of O 


These fractions contain oxygen in 
H groups absorbing at 2.95 yw (this 
band could also be due to N—-H groups), C=O 
groups absorbing at 5.86 uw, and perhaps C—O groups 
contributing to the general absorption from 8 to 
10 uw. Sample I has more C=O, O—H, C—O, aro- 
matic, and branched aliphatic groups than samples 
I] or III. Sample IIL has the least absorption by 
C=O groups, whereas sample II has the least 
branching and absorption by O-—-H groups. Sample 
I] has a stronger band at 6.21 » but weaker bands 
at 12.24 and 13.35 uw than sample III. In none of 
these are the differences in the aromatic bands very 
great 

As expected, the dark-oil fractions separated from 
the weathered asphalts show increased absorption by 
C—O groups, indicating increases in carbonyl con- 
centrations by factors of about 1.4 for sample II and 
about 2 for Land II]. There are no changes in the 
QO—H bands, and only sample I shows increased ad- 
sorption in the C—O band region. The aromaticity 
of the samples is not greatly changed, but may be 
slightly reduced by weathering. In all of the samples 
the pair of bands at 13.68 and 13.87 w have been 
greatly accentuated. As in the case of the water- 
white oils, the dark oils separated from weathered 
asphalts I] and IIT are nearly identical, although 
the dark-oil fractions from the unweathered asphalts 


are different. 
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spectra of dark-oil fractions of asphalts 


b, weathered 


III, 
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Pigure 3 Infrared spe 


3.3. Asphaltic Resins 
The asphaltic resin fractions (fig. 3) appear to be 
more aromatic than the dark oils, and to contain 
more C=-O groups. The broad absorptions neat 
3 u extend to longer wavelengths than thev do in 
the dark oil, and there is a definite absorption aut 6 pw, 
especially noticeable in sample I, suggesting the 
presence of amides or amines The strong band 
found in all samples at 9.62 uw is puzzling It occurs 
in a region where certain C—O, SQ, and Si-O 
Because the Intensity of absorp- 
these groups Is fairly large, the band at 
9.62 uw does not necessarily indicate a high concen- 
tration The band does not coincide with the band 
found in this region in the spectrum of fuller’s earth 


froups absorb 5 
tion by 


lO) gw 

Of the three samples, material | is probably most 
aromatic, and contains more OX\ venated 
than the others The band at 9.62 uw is strongest in 


samples | and III 


groups 


The fractions obtained from weathered asphalt 
show increased C=O absorption, as well as stronge! 


| hee of } ( II 
t-of-pla t on bands t veen Il and 14 mw are w 

dark-o wtion ug ting that tl isphaltic 
hly substituted aromatic rir It tot ni ! 5 t 
ilusions ba t umpt i } 

able until j uch t 4 tabl per 


ira of aspha fie sin irae Ons 
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eal: ) we 


Il 


absorption in the C Oregion. The gain in concen- 
tration of carbonyl groups measured at 5.86 uw, under 
the assumption of constant absorptis itv of carbony! 
groups, is again larger for sample | (about 2.8 times 
than for samples IT and IIT (about 1.7 and 1.8 times, 
respectively The band at 9.6 uw seems to be weaker 
in the fractions from weathered asphalt, whereas 
bands at S.S and 7.7 This is 
consistent with an oxidation of sulfoxides to form 
sulfones, for example. Once more, the near identity 
of the samples from weathered asphalts I} and Tl 
is striking. 


u become stronger 


3.4. Asphaltenes 


The spectra of the asphaltenes (fig. 4) have a 
washed-out appearance, that ts, the contrast between 
background and absorption peaks Is relatively low. 
This is often the case when rather large particles are 
observed in suspension, even though seatter is con- 
trolled by proper choice of suspension medium {7,8 
The spectra show evidence of aliphatic and aromatic 
groups and for C=O and C—O groups. However, 
no O-—-H or N—-H absorption is observed There ts 
a little difference among the spectra of unweathered 
asphaltenes and among the weathered 
specimens. Weathering appears to increase some- 
what the aromatic and carbonyl content of the 
material 


even less 
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4. Discussion and Conclusion 


Infrared analysis is most valuable when used in 
conjunction with other chemical and physical infor 
mation. The discussion of the structures found in 
the fractions ef asphalts | and IT is for the most part 
In qualitative agreement with the elemental analyses 
9) of the same or, in the case of weathered material, 
similar fractions 

All of the fractions of asphalt | have a higher ear- 
bon-to-hyvdrogen ratio than the fractions of II, except 
for the asphaltenes, and all the fractions of asphalt 
| contain more oxvgen, sulfur, and nitrogen. The 
earbon-to-hvdrogen ratio increases regularly in’ the 
order water white oils, dark oils, asphaltic resins, 
asphaltenes, and, although some oxvgen and nitrogen 
are reported im the dark oils, these elements are more 
prevalent in the asphaltic resins and asphaltenes 
The sulfur analysis, unfortunately, does not neces- 
sarily support the possible assignment of the 9.6-u 
band in the asphaltic resins to S=© groups. 

In general, weathering uppears lo reduce slightly 
the carbon-to-hyvdrogen ratio of all of the fractions 
and to inerease the oxvgen content of the dark oils, 
asphaltic resins, and asphaltenes. 

In conclusion, the information derived from the 
infrared spectra of asphalt fractions can be sum- 
marized as follows 


|. The fractions from asphalt IIL resemble those 
from asphalt Il, both before and after weathering. 

2. Water-white oils are predominantly saturated 
aliphaties, with slight aromatic content. Asphalt 
Il has less aromatic, branched aliphatic, and cvelo- 
paraffin content than I. Weathering reduces the 
aromatic, branched-chain, and evcloperaffin nature, 
particularly of asphalt | 

3. Dark oils have considerable aromatic structure 


Asphalt IL has less 


and some oxygenated groups. 
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Infrared spectra of asphaltenes separated from asphalt 1. 


, Weathered 


25-4 region, where mineral-oil suspensions were used 


aromatic, branched chain and oxygenated character 
than I. Weathering increases the carbonyl content 
and perhaps reduces the aromatic nature slightly, 
especially in asphalt I. 

$. Asphaltic resins are highly aromatic with car- 
bonvis and other oxvgen groups and possibly amide, 
S=Q, and Si-O groups; asphalt 1 containmg more 
of all of these groups. Weathering increases the 
C=O content and changes the S or Si absorption. 
Again, asphalt Lis most affected. 

5. Asphaltenes are aromatic with C=O and C—O 
eroups. Weathering produces some increases in aro- 
matic and carbonyl content. There seems to be 
little difference among materials from different as- 
phalts. 


The author is grateful to L. R. Kleinschmidt who 
provided the fractionated samples and to S. H. 
Greenfeld for his unpublished data on physical and 
chemical properties of the fractions. 
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Research Paper 2760 


Thermal Conductivity of Nitrogen from 50° to 500°C and 
1 to 1OO Atmospheres’ 


R. L. Nuttall' and D. C. Ginnings 


\ new apparatus has been constructed for measurements of the thermal conductivity 


ol gases 


ip to 500°C and 100-atmosphere pressure 
with a spacing of about 0.5 millimeter between the hotplate and coldplate. 


The parallel-plate method was used 
A capacitance 


method was used to measure the effective spacing and area of the plates under the conditions 


of the experiment 


No solid material was used between the plates. 


The effect of radiation 


Was minimized by use of polished silver parts and was accounted for by experiments with the 


conduetivityv cell evacuated 


atmospheres, 


\leasurements 
and from 50° to 500°C It is 


on nitrogen 
believed that the accuracy of the 


were made at 1, 50, and 100 


results is 


about 0.5 pereent, except at the highest gas densities 


1. Introduction 


Accurate data on the thermal conductivity of gases 
are needed for two reasons. First, thes are needed 
to cheek present theories of heat conduction in gases. 
Experimental measurements at very high tempera- 
tures and pressures are extremely difficult, so that 
theoretical means of predicting thermal conductiv- 
ities In this range are needed. Second, accurate data 
are needed on at least one gas so that engineering 
data on other gases can be obtained with relatively 
simple apparatus by a comparison method. The 
present apparatus was constructed primarily to fur- 
nish very accurate data for use as standards by 
others making thermal-conductivity measurements 
on gases. Nitrogen was chosen as the first gas to 
be measured for several More measure- 
ments have been made on nitrogen than on any other 
pure gas. It is readily available in a state of high 
purity and is entirely suitable as a standard reference 
gas for use in calibrating apparatus for relative 
measurements 


reasons. 


2. Method 


Mleasurements of thermal conductivity of gases 
ina steady state have been made by two general 
methods, radial flow and linear heat flow. 
Most measurements have been made with radial 
heat flow, either from a hot wire or between coaxial 
evliinders. In principle, these measurements are sus- 
ceptible to convection errors. Elimination of these 
errors at the high pressures requires extremely small 
dimensions, which may be difficult’ to determine. 
The linear heat flow method is free from convection 
if the heat flow ts downward, For this reason, the 
parallel-plate method was chosen for these measure- 
ments up to 100 atm. The spacing between the 
plates was made small to minimize the error caused 
by radial heat transfer. In addition, the small spac- 
ier reduces the possibility of convection in case the 
plates are not quite horizontal 


heat 
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The correction for heat transfer by radiation 
between the hotplate and the coldplate was evaluated 
by an experiment with the conductivity cell evacu- 
ated. Of course, this assumes that no appreciable 
part of the radiation is absorbed by the gas in the 
conductivity experiment. In all the experiments, it 
Was necessary to know the effective spacing between 
the plates. Frequently, solid spacers (of known 
dimensions) are used between the hotplate and cold- 
plate. No such spacers were used in this apparatus 
because it was believed that the heat transfer 
through the contact area of the spacers and plates 
might be different with gas present than with the 
cell evacuated. The radiation correction experi- 
ment would thus be partially invalidated. The 
effective spacing was measured by a capacitance 
method under the actual conditions of the experi 
ment. This method has several advantages over the 
usual method with measured spacers. First, the 
above-mentioned uncertainty in the effect of contact 
is climinated. Second, the capacitance method 
automatically accounts for any change ip dimensions. 
Third, because the direct capacitance between the 
hotplate and coldplate is measured, nonlinear heat 
flow at the circumference of the hotplate is accounted 
for. 

The thermal conductivity, & (w em7! deg C7), for 
heat flow between parallel plates is given by the 
equation 
VAN 
AtA ’ (1) 
where @ is the rate of heat (watts) flowing only by 
conduction from the hotplate through the gas to the 
coldplate, Af is the temperature difference (deg C) 
between the hotplate and coldplate, AN is the 
effective distance (em) between the two plates, and 
Ais the effective area (em*) of the hotplate. The 
factor AVA may be considered as the constant of 
the “conductivity cell” determined by the capaci- 
tance method, so that 


0.0885516Q 
('At 


(2) 





where Cis the direct capacitance (micromicrofarads 

between the hotplate and the coldplate, assuming aA 
the material between the plates to have a dielectric , a 

constant of unity corresponding to a vacuum, With | b—2 * 

the gas between the plates, a small correction must 
be made for its dielectric constant. The conduction 
equation (2) would apply, of course, to any two 
surfaces, as well as to parallel plates, provided only SS 
that the temperature and electric fields are geo- ee 
metrically similar 





ZF) 


3. Experimental Procedure 


3.1. Apparatus 


The general assembly of the thermal-conductivity | § 
cell, the pressure vessel, and the surrounding furnace | § 
is shown in figure 1. The hotplate (H), coldplate | § 
J), guard (F), and auxiliary guard (E) are shown | § 
schematically. The pressure vessel (G) is made of | } 
stainless steel and is sealed with a Monel gasket (D). | § 
This vessel is surrounded by a furnace (KK) and 
furnace “neck” (C), which are made of aluminum 
and are equipped with electric heaters in numerous 
porcelain tubes in the alumium. The furnace 
temperature was automatically controlled by means 
of a platinum resistance thermometer and bridge 
eimcuit The furnace neck was controlled relative 
to the furnace by means of a thermocouple. The 
pressure vessel extends upward out of the furnace 
region so that the electric leads can be brought out i 
of the pressure vessel in a cold region. The cooling 
coil (B) dissipates the heat from the furnace, so that 
the Lop Is cool The electric leads all No we cold 
wire) are brought up from the auniliaryv guard 
through three Inconel tubes, which serve as electro- ny pd Se ee eee * 
static shields, as deseribed later. There are 29 of |! 
these leads, which are brought out of the pressure 
vessel through a pressure seal, A. These leads go 








out radially between two “Kel-F” (polychlorotri- : 
fuoroethvlene) disks, which are pressed together for 
the seal C 
A number of difficulties were encountered before 
obtaining a successful seal at A At first, the pe a 
material Teflon (polvtetrafluoroethyvlene) was used, 0 RD 
but it was found to flow excessively at the high E 
pressures necessary for the seal. The method 
finally used for this seal was to mold Kel-F around ‘ 
the cold wires at about a ¢* G 


The vital paris of the thermal-conductivity cell 
are shown in figure 2. The hotplate (AL), coldplate L 
QO), guard (EF), and auxiliary guard (B) are all made 
of silver to minimize temperature gradients and heat 
transfer by radiation. The hotplate is made of three 
parts, silver-soldered together. The hotplate heater Q ce, A ie 1 Sy 
hs Is located between the lower parts and cOnsIsts of 





» 


about 55 ohms (at 25° () of (0.05-mm diameter O 0.5 O INCHES 
platinum wire insulated with mica. Gold leads ‘ : 

from this heater are brought out through the tem Fictre 2. Thermal-conductivity 
pering region located at HT between the upper two |, —— ata 

silver pieces of the hotplate and then to the thermal tint tance thermometers: | 

tie-down Ke The purpose of this thermal tie-down - 4 7 iad : 
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is to bring the leads as close as possible to the tem- 
perature of the silver. This type of tie-down 1s used 
extensively throughout the apparatus and consists 
of a gold terminal insulated between thin mica disks 
and held down to the silver by a nut threaded on a 
machine screw. The gold terminal is spaced around 
the screw so that it is insulated from it. 

The hotplate is supported by the guard ring at 
three points, using quartz supports (J), so that che 
bottom plane of the hotplate is in the same plane as 
the bottom of the guard ring. Because the 
efficient of thermal expension of quartz is lower than 


CO- 


that of silver, an aluminum insert (K) having a 
higher coefficient than quartz was used with the 


quartz for compensation. The quartz is held tightls 
against the hotplate by the silver screw (G) pressing 
against the aluminum insert. The guard (E) sur- 
rounds the hotplate, except where it is exposed to 
the coldplate. In the vital region near the cold- 
plate, the guard is spaced about 0.5 mm from ‘he 
hotplate. ‘The assembly of the guard and hotplate 
is held about 0.5 mm from the coldplate by three 
quartz spacers (N) between the guard and coldplace. 
In this region, the silver was highly polished to 
reduce heat transfer by radiation. Heaters (not 
shown) are insialled in both the guard and coldplate 

The temperatures of the guard and coldplate are 


measured by using the platinum resistance ther- 


mometers (D) and (Q), respectively, in conjunetion 
with a Mueller bridge. These strain-free  ther- 
mometers are not sealed. All leads (gold) from 


these two thermometers are thermally tied down at 


Cand P, respectively, so that the temperatures of 
the platinum resistance thermometers are not 
affected by heat conduction along the leads. The 


temperature difference between the guard and hot- 
plate is determined by means of a multiple-junction 
thermopile, using four No. 36 Chromel P-Alumel 
This thermopile has one 
the 


thermocouples in. series 
set of electrically insulated junctions at F on 
hotplate and the other set at C on the guard. 
Experiments indicated the desirability of providing 
another thermal tie-down zone for all electric leads 
before thes vo oul from © to the cold region For 
this purpose, the silver auxiliary guard (B) was used 
with a heater built in to provide the bulk of the heat 
flow up along the eleciric leads. All three silver 
B . (Kk ana () are held together by long 
machine serews fastened in the bottom of the sur- 
These four pieces are all 


pieces, 


rounding pressure vessel 
electrically insulated from one another 

The effective operation of the thermal-conductivits 
apparatus depended upon the proper temperature 
control of the For this pur- 
pose, five automatic thermoregulators were used to 
control the temperatures of the coldplate, guard, 


Various components 


auxiliary guard, furnace, and the furnace neck. 
These thermoregulators were actuated by thermo- 
couples or resistance thermometers (not the meas- 
uring resistance thermometers), and thev consisted 


essentially of “chopper-amplifiers” operating satur- 
able reactors in the heater circuits. The guard was 
maintained at the same temperature as that of the 


hotplate by means of the differential thermopile pre- 
viously mentioned. This temperature control was 
the most important; it maintained constancy of 
temperature to about 0.001 deg C. The auxiliary 
guard was maintained at the same temperature as 
the guard by a single differential thermocouple. As 
a result of the effectiveness of the thermoregulation 
in the thermal-conductivity apparatus, all the ob- 
servations on the thermal conductivity of nitrogen 
were made by one person. Both the electric power 
in the hotplate heater and the thermocouple readings 
were observed with a precision potentiometer. 

A capacitance bridge was used to measure the 
capacitance between the hotplate and the coldplate. 
This bridge permitted the accurate measurement of 
small direct) capacitances, even though relatively 
large capacitances exist between the “ground” and 
the plates. To eliminate, as far as possible, errors of 
measurement in capacitance, a calibrated capacitor 
having about the capacitance of the thermal-con- 
ductivity cell (about 10 wuf) was used as a reference 
and the measurements made by a_ substitution 
method. In this measurement, the same internal 
capacitance to ground was kept in the bridge circuit, 
so that the capacitance measurement approached 
very closely to the ideal comparative measurement, 
giving an accuracy comparable with the accuracy of 
the calibrated capacitor. Much care was taken to 
avoid significant capacitance between the leads from 
the hotplate and coldplate. They were effectively 
shielded from each other by bringing leads out 
through separate grounded Inconel tubes and by 
proper location of the leads that pass through the 
pressure seal, 

For the measurement of the pressure of the gas in 
the conductivity cell, calibrated Bourdon gages were 
The requirement for accuracy on these gages 
Was not great because the variation of thermal con- 
ductivity with pressure was much smaller than the 
change with temperature, 


used, 


3.2. Purity of Nitrogen 


The nitrogen used was extra-dry high-purity, 
obtained from the Linde Air Products Co... who 
stated that impurities did not exceed 100 ppm. 


Cryoscopic measurements made by G. F. Furukawa 
at the National Bureau of Standards with an adia- 
batie calorimeter showed  liquid-soluble — solid-in- 
soluble impurities to be about 10 ppm. The conduc- 
tivity cell was filled with this nitrogen, which had 
been passed through a silica-gel drver and a filter 
made of glass wool, 


3.3. Procedure 


The general procedure was to set the furnace con- 
trols to the desired temperature, fill the conductivity 
cell with gas at the desired pressure, and put in a 
chosen electric power in the hotplate. With the cold- 
plate automatically regulated to a constant tempera- 
ture and the guard temperature automatically con- 
trolled to the temperature of the hotplate, the hot- 
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plate temperature became constant. After 
capacitance between the two plates was measured, 
alternate measurements of electric power (W) and 
temperature difference (Af) between the plates were 
made until the constaney of their values indicated 
that a steady state had been reached. The power in 
the hotplate was then changed so that similar meas- 
urements were made at three or more powers. The 
gas pressure was then changed and measurements 
were repeated for each of the pressures 0.7, 50, and 
100 atm. All of these measurements were made at 
50 100’, 200°, 300 100°, and 500 

According to eq (2), the three quantities needed to 
determine the thermal-conductivity coefficient, &. 
are the rate of heat flow, Q@, the temperature differ- 
ence, Af, and the direct capacitance, C, between the 
hotplate and coldplate. Each of these quantities is 
discussed in detail below, 

In order to obtain Q@ from W, the measured electric 
power in the hotplate, it is necessary to account for 
all heat flow from the hotplate other than by conduc- 
tion through the gas. The possible paths of this heat 
flow are (1) to the surroundmg guard by conduction 
and convection through conduction through 
lead wires, and radiation, (2) to the coldplate by cou- 
vection, and (3) to the coldplate by radiation 

Preliminary experiments were made to determine 
the heat-transfer coefficient between the guard and 
hotplate. Results of these experiments showed that 
an uncertainty in conductivity of less than 0.1 per- 
cent would be caused by a temperature difference 
between the cuard and the hotplate of 0.001° C It 
was found possible to control the temperature dif- 
ference automatically to about this value 

The conductivity cell was designed to eliminate 
any convection currents. If convection does exist, 
it will cause an apparent change in the value of & 
with change in power input. At the lower gas 
densities, no such change was observed, and con- 
vection effects seemed absent. However, at the 
higher gas densities, # appeared to vary with power 
The apparent k& values were corrected by extrap- 
olating to zero power input. This apparent change 
in & with power is believed to be due to a “chimuey”’- 
tvpe convection resulting from gas flowing mto the 
space bet ween the guard and coldplate and out of the 
holes (not shown tn fig. 2) provided for electric leads 
in the top of the guard. It is expected that blocking 
of these spaces would have eliminated this effect 

The heat transfer by radiation was accounted for 
by measurements with the cell evacuated. This 
power (W,) transferred between two parallel plates 
can be expressed is 


W,—eAo(T$ 


as, 


T?), f 


where ¢ is an effective emissivity of the surfaces, <1 ts 
an effective area, o is the Stefan-Boltzmann constant, 
ana (i and T. are the absolute temperatures If the 
temperature difference, Af, is small compared to the 


average temperature, 7=—0.5 (7,4 7 then the 
transfer equation can be simplified to 
i teAoT Af BT Af i" 


the 





The constant B was evaluated in preliminary eXxperi- 
ments over the entire temperature range with the 
conductivity cell evacuated. The gas pressure in the 
cell in’ these experiments was estimated to be less 
than 107° mm of mercury. During the progress of 
the experiments with gas, the constant B was checked 
periodically at temperatures of 350° C and below to 
test for possible changes in emissivity with time 
Higher temperatures were not rechecked for reasons 
to be discussed later with temperature measurement. 

The direct capacitance ((') between the hotplate 
and coldplate was measured by a substitution method 
using a calibrated standard and a Sylvania type 125 
capacitance bridge. The precision of the  bridee 
reading was better than the certified accuracy of the 
standard capacitor. . 

The capacitance, C, in eq (2) is that measured jn 
vacuum. For measurements with gas in the cell. 
correction must be made for the dielectric constant. 
the gas. ¢€ was determined by using the 
Clausius-Mosotti equation . | 


e, of 


: = Dp, Hi 

eT" & 
where 77 is a constant independent of temperature 
for nonpolar gases, and p is the gas density. The 
value used for D was 19531077 [1] ? when p is in 
Amagat units [2]. Although the correction for the 
dielectric constant of the nitrogen may amount to 
1 or 2 percent at the highest densities, it is believed 
that the uncertainties in & due to this correction are 
less than 0.01 percent. 

Temperature measurements were made with thi 
resistance thermometers in the coldplate and the 
guard, together with the four-junction thermopile 
between the guard and the hotplate. 
calibrated by 


These resist- 
ance thermometers the NBs 
Temperature Measurements Section before assembly 
of the cell It was found, however, that when the 
system was kept under high vacuum at temperatures 
100° C, the 
changed slightly. This was probably due to con- 
tamination of the platinum by other metals in the 
Such a change took place during the deter- 
mination of radiation corrections at 500°C. Rather 
than dismantle the recalibrate the 
thermometers, it was thought that sufficient accuracy 
the 
In this way, the temperature 


were 


above thermometer characteristics 


svstem. 
apparatus to 


could be retained if two thermometers were 
compared in’ place 
difference would be known more accurately than the 
absolute temperature. This procedure is velid only 
because the evaluation of the temperature difference, 
Af, is the important factor, whereas the change in 
thermal conductivity with temperature is small. It 
is estimated that after this procedure, the 2 ‘CUracy 
of the measurement of Af was 0.002 dez C and of 
average temperature 7 was 0.01 des C 
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4. Results 


The results of the measurements are given in 
table 1. The first column gives the temperature of 
the experiment, which is the mean of the hotplate 
and coldplate temperatures. The column 
e1ves the power i put imto the hotplate (as mecas- 
ured electrically with nitrogen in the conductivits 
cell. Values of power W,, given in the third column, 
are based on experiments with the conductivity cell 
evacuated and are calculated from eq (5). These 
values are subtracted from the corresponding values 


second 


of Wo to give , as used in eq (2). Over 40 experi- 
ments were made to determine the value of B in 
this equation A value of B=4.90™ 10 w/deg in 
this equation was found to fit the results best. The 
values of B ealeulated from the individual vacuum 
experiments deviated on the average about 4  per- 
eent from the above value, ana seemed to be inde- 
pendent of temperature. Because the 
values of W) are only 3 percent of the corresponding 
values of Wy an error of 4 percent in B would cause 
2 maximum error of only about O.1) percent in 
conductivity &. 


Marini 


The values of temperature difference Af given in 
column 4 are based on measurements with the two 
resistance thermometers and the thermopile, Which 
“null indicator. Column 5 
elves the values of capacitance (© (as used im eq 4 
between the hotplate and coldplate as determined 
at the same time and under the same conditions of 


was used only as a 


the conductivity. experiment. These values have 
been corrected for the dielectric constant of the 
nitrogen 


The values of “apparent” thermal conductivity, 
i*, given in column 6, were calculated from eq (2). 
These values, when corrected to even temperatures, 
seem to be free from error due to convection in that 
they are usually independent of power, except at 
some of the highest pressures, as mentioned earlier. 
Because this effeet is believed to be due to a chimne\ - 
type convection, it should diminish at lower power 
inputs. Therefore, the values of 4* were extrapo- 
lated to zero power to give “true’’ conductivities, & 
In the worst case, at 100-atm pressure and 50° (©, 
the extrapolation is over a large range of values of 
k*. However, it is believed that this extrapolation 
should introduce only a relatively small error. The 
results from the three experiments at different powers 
were fitted by the method of least squares to a 
linear equation in power. The maximum deviation 
of the observed values of /* in this linear equation 
was only 0.28 percent, so that it is believed that the 
linear extrapolation was adequate. 

At low pressures the effect of temperature discon- 
tinuitv near the walls may become significant. This 
has been taken into account for the 0.7-atm data by 
ise of the equation 
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The “temperature-jump distance’, g, 1s 
the “accommodation coefficient’, a, by the relation 
116): 

2—4a,, in k 


DOP 


T’ the absolute tem- 


Ccon- 


in which # is the gas constant, 
perature, & the thermal conductivity, Cy the 
stant-volume specific heat, y the specific-heat ratio, 
and P the pressure of the gas. Values of a for nitro- 
gen on silver vary from about 0.8 at room tempera- 
ture to 0.4 at SOO° C [6] The value of 0.5 was 
chosen as the most probable for the temperature 
range of these experiments 

Values of observed conductivity 4 
even temperatures, are listed in column 2 of table 2 
In order to facilitate interpolation of the results be- 
tween the observed points, empirical equations were 


corrected to 


derived to give conductivity as a funetion of tem- 
perature at each of the three pressures. These 
equations are as follows 
q 
r=] 059 ( (4 h ) pressure, O.¢ atm 4 
t 
ai: 2? 714-0.005S8971 pressure, 00 atm), 10) 
Lh = S.1OL + 2.9317 © 10+ 9.0761 < 10° 
S.SS1S 107" pressure, LOO atm 1] 
In these equations, i is in watts em deg", f is in 
deg (, 9 Is: viscosity poise), €y is heat capacity 
} ¢ deg at constant volume, and / is the gas 


constant 8.517/28.016 a deg The values of 


viscosity and heat capacity used in eq (9) were ob- 


Tarte 2. Thermal conductir 
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NBS Circular 564 (Nov 


| { ’ D ra) 


' Which 
tabulates heat-transport properties for a number of 


tained from 


The constants of these Were 
obtained by the method of least squares. Values 
of & given in column 3 of table 2 are calculated from 
these equations. Column 4 gives the deviations 
of the observed results from these calculated Values 
The results are viven in figure 3, where the 
circles represent observed values and the solid lines 


equations 


CAaSCS. 


also 


the values from the above equations 

The empirical eq (10 and (11) fit 
Values at 50- and 100-atm = pressure, respectively 
almost as well as would be expected from the pre- 
cision of the data The 
1O)(which is linear in?) from the observed condue- 
tivities is only 0.1 percent, with a maximum deyi- 
ation of 0.3 percent, The deviations from eq (11 
are also quite small except below 200°C wher 
there was increased difficulty, both experimentally 
and in fitting the observed points with an equation, 
The deviations of eq +t from the observed Values 
seem to be larger than can be attributed to experl- 


the observed 


average deviation ol eq 


mental errors. Equation (9) is the same semi- 
empirical equatl rn proposed bry Kucken ot, except 
for the constant factor A faetor of 1.059 was 


found to fit the present data, as compared to a factor 
of | originally proposed Iyy Kucken Nore recent 
modifications of Eueken’s equation by other in- 
vestigators eenerally have predicted i value slightly 
vreater than 1. An alternative equation for 0.7 
atm, which fits well up to 400°C, and which may bi 


easier for use in numerteal interpolation in this 
region, Is 
/ 2 495 10) OH 306 1Q-f OOD 10) » 

. P=100 atm 

*k=3.16 + 2.9317 X 
o +9.0761 X 107®¢? 
C § 8.8318 x 107? t3 - 
SO atm 
> kk 7ia+ 589 + e 

2 

x 45 
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where & is in w em™! deg™', and ¢ is in deg C. Al-- 
though all of these equations are suitable for inter 
polation at their respective pressures, they may not 
be valid for extrapolation bevond the experimental 
range. Equation (9) is probably the most suitable 


for extrapolation to higher temperatures. 


5. Discussion of Results 


5.1. Accuracy 


estimate the acecuracys of the final 
to consider the effect of each 


In order to 


esults, it 
fthe terms in the conductivity equation 


Is Hhecessars 


0.08855] it) 


(‘Al 
The consiant O.OSS5516 is a conversion factor to 
allow for the svstem of units. It has negligible 
neertainty in its use here. Some possible sources 


capacitance between electric leads and the relatively 
large capacitance to ground (the guard is grounded ) 
of the hotplate and coldplates. The electric leads 
were carefully shielded to avoid significant error 
from this source. The effect of the ground capac- 
itance is largely eliminated by the design of the 
bridge used in its measurement. Use of a substitu- 
tion method that compares the conductivity-cell 
capacitance with that of a standard capacitor 
minimizes a number of possible errors. By using 
this method, it is believed that this comparison is 
accurate to better than 0.1 percent. The absolute 
value of the standard capacitor is certified by the 
NBs (with a high probability) to only 0.3 percent. 

fn the experiments wiih nitrogen, the effect of 
dielectric constant of the must be con- 
sidered. It estimated that uncertainty in the 
value of the dielectric constant of nitrogen is less 
than 0.01 pereent, giving the same uncertainty in 
Cand k The total unceriainty in the capacitance 
measurement is believed be about 0.3 percent 
This is the largest uncertainty in most of the meas- 
and it is based on 1 chance in 10 that 


vas 


Is 


urements, 


of errorjin the measured capacitance, C, are the | the error will be larger than 0.3 percent. 
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The factor Q in eq (2) is evaluated from the two 
powers Woand W,. The electrical measurement of 
W is believed to have negligible error. The ulh- 
certainty in the value of W, determined from vacuum 
experiments is estimated to be less than 0.1 percent, 
except at the highest temperatures, where it| may 
be as large as 0.1 to 0.2 percent It is possible that 
the vacuum experiments do not account completely 
for all heat flow (other than Q) from the hotplate 
during the gas experiments. This may be due to 
imperfect matching of the guard temperature to the 
hotplate temperature. It is that the 
error from this source is less than 0.1 percent 

The measurement of At, which has been discussed 
earlier, has an estimated uncertainty of 0.002 deg 
(. This introduces uncertainties in the values of 4 
of from 0.02 to 0.2 percent, depending on the mag- 
nitude of Af. The uncertainty in the temperature 
scale, estimated to be about 0.0] deg = believed 
to introduce negligible error in & 

In estimating the over-all accuracy of the results 
the authors believe that with most of the results, 
there is only 1 chance in 10 that the true conductivity 
values will deviate more than 0.5 percent from the 
observed values fiven in table 2 In the case of a 
few values at the highest believed 
that this tolerance might 
as | pereent 


estimated 


pressures, it is 


he increased bry iis much 


5.2. Comparison With Others 


\ comparison of the results of this mnvestigation 
with the observed conductivits values from some 
other recent given in figure 4 No 
attempt was made to include in the figure all the 
available data at low temperatures and low pressures, 
although all known high-pressure data are vive 
The data are plotted as deviations from eq (9), (10 
and (11), which represent the NBS data for the three 
pressure ranges At low | 


researches Is 


pressures (approximately | 
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atm), Frank [4] and Schottky [5] made measurements 
up to 500° C, Rothman and Bromley [6] up to 809 
C’, Stops [7] up to 1,020° C, Michels and Botzen fg 
up to 75° C, and Vargaftik and Oleschuk [9] up to 
541° C. Keves and Sandell [10] have also mad 
measurements up to 100° © at various pressures 
Their results are not shown here because they beliey, 
their conductivity values to be too low over most of 
the temperature range. Later observed values of 
Keves [11] upto 150° C are shown. Further measure- 
ments by Keves [12] have been made up to about 
350° C, but no observed values of conductivity hay: 
been published At higher pressures, Michels and 
Botzen [S] made measurements up to 2,500 atm and 
Za” 6S; Stolvaroy, [pat’ev, and Teodorovich 
to 500 atm and 300° C, Vargaftik [4] up to 90 atm 
and 62° C, and Lenoir and Comings [15] up to 200 
atmat41°C. Asin the low-pressure measurements 
the results of eves and Sandell are not 


13] up 


show 1 


although later results of Keves ure given up to LOO 
atm and 50° CC 
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Apatites Deficient in Divalent Cations 


Aaron S. Posner' and Alvin Perloff 


Phe preparation and results of a study of stoichiometric and divalent cation-deficient 
defect) lead and ealeium hydroxvapatites are reported \-ray diffraction and index-of- 
refraction data show that it is possible to have divalent cations missing from ‘‘ecolumnar”’ 
positions the hydroxyapatite structure, with possible hydrogen bonding in proper pro- 
portion for electrical neutrality \ general formula for apatites is suggested to be 
1) oi PO) OH)... where D re presents the divalent cation (Pb*?, Ca Mg*?, ete and 
\I the monovalent cation, such as Na* and H*, and the upper limit of 2 is probably 2. It 

sugeeste d that the apatite portiol of bone at ad tooth mineral May be some form of defect 
tporuatite 


1. Introduction 


{ number of finely divided calcium phosphates 
vive the X-ray diffraction pattern of hydroxy apatite, 
Cay (PO,),(OH but depart from the ideal Ca/P 
ratio of this mineral The formulas in this papet 
are given for unit-cell contents.) Examples are to be 
found in the inorganic portions of some bone and 
teeth, the mineral phosphorites, and the chemically 
precipitated synthetie apatites [12,3] There are 
two conflicting explanations for a low Ca/P ratio 
n certain high-surface-area hvdroxvapatites. The 
first states that an excess of phosphate ions ts ad- 
sorbed on the surface of tiny erystals of stoiehio- 
metric hvdroxvapatite o4 The second explana- 
tion postulates a statistical absence of calcium ions 
from structural throughout the ervstal 
56 In this ease, electrical neutrality is probably 
obtained by the inclusion of 
vdrow 1 bonds between OXVeens ol orthophosphate 


positions 


protons present as 
groups ordinarily bonded to the missing ealeium 


data to the latter 


This paper presents support 
theory and enlarges this view to include apatites con- 
taining other divalent cations, such as lead ions 
These compounds with Imperfect C'n P or Pb P 


atios, are referred to as defect apatites In this 
this ervstallographie term implies that certain 
space-group positions normally containing divalent 


cations (Ca, Pb) are statistically unoccupied 


isnle 


2. Adsorption Theory 


\lany workers have said that the large surface area 
of the chemically precipitated, “apparently 
stoichiometric” hvdroxvapatites can accommodate 
CXCECSS phosphate ions neutralized yy some adsorbed 
cation [3.4 If this were true, it 
assume that the surface exchange of phosphorus com- 
pared to calerum would be higher than the over-all 
Ca P ratio Using radioisotopes, Weikel and Neu- 
man {7! incieated that this was not true The differ- 
ences in chemical analysis of apatites studied could 
hot be correlated to the surface composition. Ac- 
the low Ca’ P apatite did not 


honh- 


is reasonable to 


cording to this work. 


show more phosphate surface exchange than the 
normal Ca/P apatite. 

A few simple calculations can shed serious doubt 
upon the surface-adsorption theory. In the 
of any apatite, to lower the divalent cation to phos- 
phorus molar ratio from 10/6 to 9/6, there must be an 
adsorption of 0.67 mole of phosphate per mole of 
apatite [S}. Considering each tetrahedral PO, group 
as a sphere of radius 2.96 A, the effective area occu- 
pied on the apatite surface by each PO, group would 
be 27.53 A?. 

Because 1 g¢ of apatite is about 10 mole, then 
0.67 > 107° mole of surface phosphate ion would be 
needed for 1 g¢ of this compound. Therefore, 1 g of 
apatite would require a surface area of 111.1 m* to 

the above amount of phosphate 
97.53 X10 6.02 “10° 0.67 10 111.1 m*/g). 
If the surface diminished below this value, the phos- 
phate could not be adsorbed in great enough quan- 
tity to lower the CaP to 9/6 from 10/6. Both defect 
calcium hvdroxvapatite and defect lead hvdroxy- 
apatite have been prepared with large enough ery stal 
sizes to insure a maximum surface area far too small 
to account for the lowering of the metal to PO, molar 
ratio by adsorption of excess phosphate. This dis- 
cussion has been oversimplified and does not take 
into account the eryvstal surface needed to accommo- 
date any cations that may with the 
adsorbed phosphate ions. 


Case 


accommodate 


be associated 


3. Preparations and Experimental Methods 


A series of calcium and lead apatites, with both 
ideal and defect divalent cation content, were pre- 
pared. The calcium compounds prepared ranged in 
(a P weight ratio from the theoretical ratio of 2.15 
down to 1.82. The two samples of lead apatite pre- 
pared had, respectively, the theoretical Pb/P ratio 
of 11.10 and the “nonstoichiometric” ratio of 10.03, 


3.1. Calcium Apatite 


The low Ca/P (1.82) (sample 6, table 1) apatite 
was prepared according to Arnold [9]. A liter of 0.1- 
molar CaCl, was mixed with about 7 liters of 0.1- 
molar Na,HPO, and stirred for 2 days at room tem- 
perature. The solid was filtered, washed free of Cl- 
and Nat, dried at 110°C for 24 hr, and analyzed 
chemically. 
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The apatitic caleium phosphate with Ca/P ratio 
of 1.94 (sample 5) was prepared by the method of 
Dallemagne [10]. A saturated solution of Ca(OH), 
was added to a solution consisting of 4 ml of S85- 
percent H,PO, in 4 liters of water until the pink end 
pomt of phenolphthalein was visible. Filtermg and 
drvine (110° C, 24 hr vielded a defect apatite whose 
molar ratio was 9 Ca to 6 P. 

The other calcium compounds, which had Ca P 
ratios 2.08, 2.11, and 2.15 (samples 4, 3, and 1, re- 
spectively), were all a result of the Wallaevs {11 
procedure for the preparation of synthetic apatites 
The method is similar to the Dallemagne 
method, but after neutralization of the acid solution 
by Ca(OH),, the H,PO, solution is boiled and neu- 
tralized at the boiling temperature by the 
method. If the neutralization point is reached, the 
product is an apatite of theoretically correct Ca P 
ratio (2.15). By falling short of the neutralization 
point at boiling temperatures, the lower Ca P values 
of 2.11 and 2.08 were obtained. 

{nother method was used to obtain stoichiometric 
calcium hvdroxvapatite (sample 2) [12]. Pure crys- 
talline CaHPO, was prepared and hydrolyzed to 
hvdroxvapatite with distilled water im a platinum- 
lined hydrothermal bomb. This method vielded 
larger ervstals than the Wallaevs preparation 


above 


same 


3.2. Lead Apatites 


Pure lead apatite [Pb (PO,),(OH 
was made by placing S ml of 0.5-17 lead acetate 
solution ina platinum-lined hvdrotherma! bomb and 
slowly pipetting 4 ml of 0.6-\7 (NEH,).HPO, solution 
onto the first laver. Care was taken to maintain a 
The bomb Wits sented anid held 
which time the svstem was 
The resulting well-ervs 
filtration, afte 


sample 7 


two-laver system. 
at 300° C for 7 davs, at 
quenched rapidly in wate 
tallized solid was collected by 
it was washed in water, then in acetone 

\ commercial lead phosphate (Fischer Scientific 
1.5), whieh had a molar Pb-to P ratio of 8.9 to 6 
and gave the A-rav pattern of lead apatite, served 
us2 sample of defect lead hvdroxvapatite sample ‘ 
According to the manufacturer, it was prepared by 
mixing solutions of tribasic sodium phosphate and 
lead nitrate and decanting the supernatant liquid 
The precipitate was washed with distilled water by 
decantation until the test for nitrate and sodium was 
hevuative The material was then filtered and dried 
110° C, 24 hr 

In all the anmalvses: (a) caletum 
volumetrically by titration of the 
standard potassium permanganate 
ly phosphorus Wis determined colorimetrically by 
the heteropoly blue method [14], and (e 
determined gravimetrically as PbSO, 


3.3. Bomb Growth 


chemically precipitated apatites are usually 
divided materials vielding X-ray diffraction 
In two cases the 


which 


was determined 
oxalate with a 
solution {13}, 


lead Wiis 


The 
finely 
patterns with broadened maxima 
apatite crystal size was increased by hydrothermal 
bomb vrowth 15, 1S] 





Two grams of low Ca/P calcium apatite sample 5° 
were introduced into the platinum-lined chamber o* 
the bomb and covered with 4 ml of distilled water 
The bomb was sealed with a platinum gasket and 


inserted into a furnace at 300° C' for 7 davs. The 
resulting solid was filtered, washed, and dried, and 
designated as sample ja 

In the ease of the lead sample, 0.52 g of sample s 
and 10 ml of distilled water constituted the bomb 
charge. The bomb was held at 300° C for S days: th, 
sumple was removed by filtration and designated as 
sample Sal 

In both CUuses the solid phiase did hot change in 
chemical analvsis after the bomb treatment In 
the case of the ecaleium apatite erowth the bomb 
water showed only trace quantities of caleium and 
phosphate lions 


3.4. Physical Measurements 


The average indices of refraction of the prepara- 
tions were obtained by using oil-immersion  tech- 
niques with the polarizing microscope. X-ray 
der diffraction studies were carried out on 
tional-counter X-ray spectrogoniometer 
nickel-filtered copper K-radiation, taking 
precautions in preparing the sample to avoid erystal 


pow- 
a propor- 
using 


special 


orrentation 

In the detailed N-rav study of the line intensities 
of the two lead apatite samples, a digital recorder was 
attached directly to the X-ray scaler and set to print 
out the integrated intensity every LO see. Beecausi 
the traveled at) 0.2° 20) pet 
minute, a LO-see print-out corresponded to an ints 
vration over a range of 0.033° 20. The intensity of 
each diffraction maximum was caleulated from. thy 
digital Information appearing on the ndding-ma- 
chine tape print-out svstem 


spectrogoniometer 


4. Results 


thie results obtatmed for thi 


Table 1 
ealerum upatites 
lt is possible to caleulate by thy method of lorentz 


prese! i 


and Lorenz 16] the theoretieal change In average 
index of refraction as ecaletum ions are removed 
statistically from = structural positions. Figure | 


eompares the variation in average index of refraction 
in calcium apatite as the Ca P is varied by missing 
calcium ions (data for sample 2 omitted 

Table 2 presents the results of the studies on the 
lend npatites Included in the table is a Comparison 
of the observed and caleulated index of refraction for 


both the stoichiometrie and defect materials The 
Gladstone-Dale [17] method was used for thes 
calculations 

The X-ray diffraction patterns of the caleim 
apatite samples were all very similar \s reported 
earlier [S|], there is a shght tendenev for the unit- 
cell volume to shrink as the Ca P ratio falls How- 


ever, the relative intensities of the maxima do not 
change appreciably with changing calcium content 
On the other hand, there are detectable differences 
in the relative intensities and spacings of the N-ra\ 
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Data 


NS ee Ee en Se 5. X-ray Study 
5.1. Stoichiometric Lead Apatite 


The powder pattern of lead apatite (table 4) was 


Average ind x 
tt found to be in good agreement with earlier reports 

(18, 19]. A struetural analysis of the data was per- 

ness | ws pen formed to see if the fluorapatite parameters were 

i 45 £0 0 Do. suitable for lead apatite. Because so many of the 

. 0.00 ID X-ray maxima of this complicated structure repre- 
er — sent the superimposition of many independent re- 

. ‘ M D flections, the X-ray analyses in this work were carried 


out by using the relative-intensity values rather than 
the structure factors. 

The space group of apatite is C%,(P6;/m). The 
atomic parameters assumed were taken from 5. 
Naray-Szabo [20] except for the z-axis parameter of 
94 910 P, which was suggested by C. A. Beevers as an im- 
proved value (private communication). 





- | LR » Liomire parameters of the asummetrir init of the 
nal ee ol tuo apatite applied to lead h ydroxya patil 
a f 
. 
Numb P meters in fraction of axis lengt! 
\l | 
. 
. 
>} j i4°8 Oo. M7 oO. O00 
Ph 4 om ooo 20) 
OH » wn mn 2) 
o Oy ‘ 33 un) 2) 
0 ) Wi7 7D 
oO 12 e°2°3 oo) (ans 
| ‘ cH) eu) oO) 
| Na Szat 2 





Using the parameters in table 3, the structure 
of refraction with CaP | faetors, and in turn the theoretical intensities for the 
observed X-ray maxima, were calculated. Compar- 
ing the observed and the calculated intensities of the 
first eight reflections (table 4), one arrives at a 
reliability factor of 21.8 percent. This is good agree- 
ment for jntensity (as opposed to structure factor) 
comparisons, and the atomic parameters of fluor- 
apatite are acceptable for this powder diffraction 
study of lead apatite. The employment of the 
multiple-index, higher-order reflections would have 
added little change to the reliability factor 
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5.2. Defect Lead Apatite 


The four Pbt (columnar) atoms are stacked in 
two columns in the cell interior on the threefold rota- 
tion axes parallel to the c-axis. The six Pbu 
(hexagonal-screw) atoms are found associated with 


the hvdroxvl ions. W hich are located Ol the hexagonal- 
screw axes parallel to the c-aNnis Figure Za diagrams 
a projection of the atoms of lead apatite upon the 


basal (001 plane Figure 3 is the same projection 


with only the lead atoms shown. 
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The molar Pb-to-P ratio in the sample studied was 
8.9 to 6, but for simplicity of calculation, it) was 
assumed to be 9 to 6, with one lead ion missing per 
unit cell. This cation may be absent in three dif- 
ferent Ways: (2) missing statistically from all possible 
cation positions; (b) L of the 4 Pb rt atoms may be 
missing statistically; or (¢) Ll of the 6 Phat 
may be missing statistically. There is also the possi- 
bilitv that a specific cation is missing in each unit 
cell, but the lack of X-ray evidence of a superlattice 
made this improbable, and the possibility. was ig- 
nored Because the above parameters were found 
acceptable for lead apatite, the structure factors for 
reflections out to sin O@rA—O.4ATL (A—1.5405) were 
calculated for the above possible arrangements and 
for ideal lead apatite. 

Table D lists the observed and cnieulated 
ties for the first 28 reflections of the defect apatite, 
As can be seen, the over-all reliability factors for the 
This is not 


atoms 


itensi- 


four possibilities are all reasonable 
SUPPPrising because the majority of the reflections are 
litthe changed by a missing lead ion. Table 6 sum- 
marizes the difference between observed and caleu- 
lated tensits 
alfected reflections 
b). in whieh one of the four Pb t atoms is MISSIN 
tuble 6, last column), is the best mateh 


values of these eight most strongly 


There is no question that case 
per unit cell 


6. Discussion 


6.1. Lead Apatite 


The X-ray data obtained from defeet lead apatite 
indicates that the cation is missing statistically from 
the Phi fie. 3 Thus, im the 
case of a “Pb P apatite, | of 4 cations located 
on the trigonal axes (i. e., Pbt) is missing 
cell, with the defeet statistically distributed through 
the } positions . 


structural position 


per unit 


index-of-refraction data for 
The observed drop ih 


This corroborates the 
lead upatite given in table 2 
index of refraction 


atom missing from the structure 


Is as expected, with one lead 
lndex-of-refraction 
data alone Cal merely Suggest a missing Cation pel 
unit cell The X-ray data identifies the exact nature 
of the defeet structure. Reeall that the problem was 
simplified by assuming one Pb ion missing per unit 
cell Actually, the chemical analysis of the 

storchiometrie apatite indicated that there was a 
defeet of 1.1 lead atoms per unit cell This empha- 
sizes the point that statistical absences of cations in 
such a structure discrete numbers 
This is well illustrated by the calcium apatite series. 


hot- 


need not be In} 


6.2. Calcium Apatite 
The index-of-refraction data (fie. 1) show that the 
calcium. series parallels the lead apatite case As 
more cations are removed, the observed index drops 
It is not surprising that the X-ray patterns of stoich- 
the nonstorhiometric calcium 
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\ caleulation of the effeet of 
cell indicated no measur- 
number of lead as 


apatites are so simula 
a missing Ca ton per unit 
able change. The high atomie 
compared to caleium results in a greater change in 
the X-ray intensities when | in 10 atoms is missing. 

It is fortunate that the missing calcium atoms do 
not affeet the ervstalline unit-cell volume appreci- 


observed X-ray powder data on defect lead apatite with the calculated intensities for the four possible cases 
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ably. As the number of calcium atoms per unit cell 
declines, it is inevitable that the index of refraction 
falls, because the unit-cell volume ‘is reduced only 
slightly. The lack of exact match between caleu- 
lated and observed indices of refraction in the cal- 
cium apatite series is probably due to the inexact 
nature of the molecular refractivities assumed for the 
l_orentz-Lorenz calculation. It is important to note 
the parallel nature of the drop in index of the ob- 
served and theoretical apatites as the caleium con- 
tent of the unit cell is reduced. It seems reasonable 
that a defect apatite series exists in the case of cal- 
cium with cations missing from columnar positions 
as in the lead compounds. Here, too, the calcium 
ions are probably absent statistically from the four 
unit-cell positions on the trigonal axes and not in any 
integral number per unit cell. 

An earlier study [8] showed that the observed 
density of a caletum-defect hydroxyapatite (with 
Ca P=9/6) was 2.89, as opposed to 3.13 for pure 
hydroxyapatite. This is the value predicted by 
assuming the missing calcium mechanism. It serves 
to corroborate the index-of-refraction data pre- 
sented above 
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6.3 Hydrogen Bonds 


Electrical neutrality in these low-cation apatites 
can be achieved by the inclusion of hydrogen ions 
as hvdrogen bonds in the proper proportion It is 
not reasonable to place two protons in a lattice posi- 
tion previously occupied by a lead or calcium ion 
The hvdrogen probably enters into the strueture as a 
hvdrogen bond between the OXVECHs of adjacent 
orthophosphate tetrahedra. It is not possible to 
prove this hypothesis with the present N-rav data 
Future single-crystal studies with X-ray and neutron 
diffraction will be needed to further establish this 
suggested mechanism of achieving neutrality. 

One can envision 
hybrids 
orthophosphate structures. In the caletum = series 
the strueture of hydroxyapatite is quite different 
from CaHPO, or CaHPO,2H.O [21,22], and it is 
unlikely that a continuous solid-solution series exists 
between apatite and either of the latter compounds 
There must be a lower limit to the amount of missing 
cation in the nonstoichiometric series for, as the 
Ca P atomic ratio approaches unity, there will be 
alterations in. the crvstal structure In the svn- 
thetic calcium apatite svstem, the lowest Ca P ratio 
9.25.24]. Whether this ts 
the limit or not is unknown at this writing. All 
attempts to precipitate apatites with Ca P ratios 
lower than 8 6 by the authors and others [9.24] have 
resulted ino the formation of an acid phosphate 


pl 


these apatites with missing 


cations as between the acid and basie 


reported to date is S6 


ise. 


6.4. Evidence for Hydrogen Bonds 


\n interesting series of experiments [25, 26] has 
viven some experimental evidence for the presence of 
hvdrogen in defeet calcium apatite. The caletum 
uptake from saturated Ca(OH), solution by finely 
divided “low-caleitum”™ apatites was observed by 
Chemical analysis of the solid and conduetivits 
measurements in the solution. Higher temperatures 
resulted in an increased calcium uptake, and the 
amount of ealetum chemically bound (i. e., Ca not 
removed by washing) alwavs brought the Ca P of 
the solid to 10.6, the ratio for pure apatite In all 
cases the conduetivity of the Ca(OH), solution 
dropped further than the Ca loss would require 
It is feasible that the hvdrogen leaving the apatite 
when calcium returns to the empty cation position 
either by diffusion or reerystallization) will form 
poorly ionized H.O, thus altering the conductivity 
expected by only considering the Ca loss. Neuman 
[27] confirms this by showing a lowering of the pH 
in solution when caleium is added to “basic calcium 
phosphates” (another term for low-calcium apatites 
He reports this to be due to the release of H,O* as 
the calcium ions return to the defect sites (see section 


6.7 
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6.5. Nomenclature 


In light of the probable hydrogen bonding it is not 
completely correct to refer to these compounds as 
nonstoichiometric hydroxvapatites. The more pre- 
cise term defect apatite is thus emploved., 

\s in other complicated series of compounds, the 
hvdroxvapatite “family” is best characterized by a 
general formula. This might be D,,_-M.,(PO,),(OH),. 
where D> represents the divalent cation (Catt 
Pott, Sr and M the monovalent cation (Nat. 
Hi*. ete In the ease of the defect lead apatite 
given above, the formula becomes Pb..,H...(PO,),- 
OH).. A sodium substituted apatite would have 
the general formula Dy) -Na-H,(PO,),(OH)., where 
Na can occupy the divalent cation positions in the 
lattice with hvdrogen bonding effecting electrical 
neutrality This paper does not deal with possible 
isomorphous substitutions for phosphate and 
hvdroxvl ions 


6.6. ““Octocalcium Phosphate’ 


Arnold [9] suggested a mechanism to explain the 
calcium-defieient hydroxvapatites that he prepared. 
He postulated that as many as 2 of the 4 columnar 
calcium ions per unit cell could be missing to form 
octocalcium phosphate, an apatite with a Ca P ratio 
of 86. To achieve this stoichiometry and maintain 
the structural sVinmetry of apatite, Arnold proposed 
that all the hvdroxyvl groups are Missing, whereas 
the calciums (missing at the same columnar c-level 
are replaced by water, aud that hydrogen ions are in 
some Way associated with the columnar calerums that 
are present. His theory allows for a series ranging 
between this low-calcium apatite end pure apatite 
He suggests that octocaletum phosphate is an In- 
finite two-dimensional complex, Ca,H¢PO,).3H.O, 
with sheets held together by hvdrogen bonds be- 
tween the oxvgens in the POs groups and the oxygens 
of the water molecules. As more columnar ealerum 
is. taken up, Arnold feels that hydrogen would be 
randomly climinated, resulting in ionie bonding 
between the sheets 

Even though Arnold's viewpoint superficially re- 
sembles the theorv of defect apatites presented in this 
paper, it is actually quite different. He recognized 
the need for an explanation of low-caleium apatites, 
but his mechanism must be reyectod ona few simple 
erounds First, it unlikely that the apa- 
tite structure will hold together with the hydroxy! 
ions missing. Granting that such a structure were 
possible, the removal of the OH ions from = their 
special positions would result in measurable changes 
in the relative intensity of many of the X-ray dif- 
fraction maxima. In addition, the mechanism of 
water substitution for columnar caletum between 
the “sheets” would certainly result) in’ expanded 
unit-cell constants. Such changes in the apatite 
X-ray patterns of these low-calcium materials have 
not been demonstrated. Arnold [24] and Engstrom 
125] reported certain low-angle, “nonapatite’? X-ray 
maxima in the study of octocalcium phosphate that 
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be fully explained. Segerman (private 


could not 
communication) observed similar low-angle X-ray 


reflections in studying dental calculus. He suggested 
they may be due to an ordered loss of calcium ions 
in the columnar position, resulting in an apatite 
superlattice. If this can be clearly demonstrated, 
it will add further proof to the existence of caleium- 
deficient hvdroxyapatites. In addition, such a 
mechanism may provide the structural link between 
the hydrogen-containing orthophosphates and the 
basic orthophosphates. 

Besides the lack of X-ray diffraction corroboration, 
Arnold’s theory does not correctly predict the ob- 
served index-of-refraction data. If missing caleitum 
ions are supplanted by water, the index of refraction 
of low-caleium apatites will not drop as observed 
fig. 1 Thus, Arnold’s hypothesis, although pro- 
vocaltive, Is felt to be incorrect. 


6.7. Calcified Tissue 


The major phase ol the iInore@anie portion of bone 
and teeth is generally accepted to be a low Ca/P 
form of hydroxyapatite [1]. Neuman and Neuman 
I], in a verv thorough review article on the mineral 
portion of bone, offer a concept to explain this lack 
of stoichiometry in bone apatite. Adopting Arnold’s 
view that water can replace caleitum in the apatite 
structure see section 6.6). they postulate missing 
ealclum ions supplanted by hydronium ions (H,O 
on the surface of the tiny apatite ervstals in bone. 
Because the surface in bone-apatite ervstals repre- 
sents such a large proportion of the solid, the aver- 
age Ca P of the material is thereby lowered. Like 
the Arnold theory, this view of nonstoichiometry in 
bone apatite is open to question on the basis of index 
of refraction and X-ray diffraction results 


0.2 


section 


The theory of Neuman and Neuman depends 
entirely upon a small ervstal size in the same way 
phosphate-adsorption theor, discussed in 
2. In contradiction to this, it is possible 
to prepare lead defect and calcium defect apatites 
with low specific areas (section 3.3 According to 
Neuman and Neuman, the particle size of a bone- 
apatite ervstal will dictate the Ca/P stoichiometry. 
In the results reported above (section 3.4), there is 
no correlation between crystal size (from X-ray 
diffraction line broadening) and Ca/P ratio. 

Surface-exchange experiments on a calerum defect 
apatite (Ca/P--9/6) showed that about 5 percent 
of the total calcium was found on the erystal surface 
27 | Because 10 percent of the total calcium ions 
were postulated to be missing from surface positions, 
15 percent of all the caletum positions in the crystal 
would have to be on the surface. Such a crystal 
would require a specific surface of about 200 m?’, or 
at least three times the value reported for the above 
defect. apatite with 5 of exchangeable 
ealcrum [27]. 

Bone-apatite ervstals range in size from about 
25 Ato about 1.500 A. 
about 200 to 300 A [28] 


as the 


section 


percent 


The surface-substitution 


with an average diameter of 


mechanism could only explain the calcium defects 
in the 25- to 100-A region because larger crystals 
cannot furnish the surface needed for this postulated 
effect. It seems improbable that the defects should 
exist only on the surface for smaller crystals and be 
distributed throughout the material for larger 
crystals. 

It is probable that the apatites in bone and tooth 
mineral are isostructural with the synthetic, low- 
calcium apatites described in this paper. It is 
further suggested that these biological apatites can 
take up calcium (and other ions that can isomor- 
phously substitute for calcium) from the body fluids 
and release hydrogen ions previously found as 
hvdrogen bonds in the solid. The differences in 
stoichiometry found in various specimens of bone 
and tooth salt may be related to changes effected 
by exposure to serum containing ions (Ca**, Mg**, 
Seve Snt*, ete.) that react with defect apatites. 


7. Summary 


Methods of preparation of both stoichiometric and 
defect lead and calcium hydroxvapatites are discussed. 

X-ray diffraction data are given to verify the 
use of the fluorapatite parameters for lead apatite in 
structure-factor calculations. 

The existence of defect lead hydroxyapatite and 
calcium hydroxyapatite is supported with X-ray 
diffraction and index-of-refraction data. An analysis 
of the X-ray data for cation-defective lead apatite 
indicates the lead ions to be missing statistically from 
the “‘columnar’’ positions. 

With divalent cations missing from positions in the 
unit cell, the electrical neutrality of apatite is 
postulated to be obtained by hydrogen bonding 
between oxygens of the orthophosphates surrounding 
the defects. 

A general formula for apatites of this type is 
suggested to be Dio ,M>,(PO,),(OH)., where D is 
the divalent cation (Pb*?, Cat?, Mgt?, ete.) and M 
is the monovalent cation (H* and Nat). The upper 
limit of « is unknown but present results suggest that 
this value is about 2. 

It is quite probable that the apatite portion of 
bone and tooth mineral is identical to the synthetic, 
defect apatites deseribed herein. 
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